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ABSTRACT 

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 
spectrograph (/Vfib C r=2400, 380< A < 1300nm), offering unique opportunities in survey astronomy. Following 
a successful external design review the instrument is now under construction with first light predicted in late 
2017. Here we summarize the science case for this unique instrument in terms of provisional plans for a Subaru 
Strategic Program of ~300 nights. We describe plans to constrain the nature of dark energy via a survey of 
emission line galaxies spanning a comoving volume of 9.3ft.~ 3 Gpc 3 in the redshift range 0.8 < z < 2.4. 
In each of 6 independent redshift bins, the cosmological distances will be measured to 3% precision via the 
baryonic acoustic oscillation scale and redshift-space distortion measures will be used to constrain structure 
growth to 6% precision. As the near-field cosmology program, radial velocities and chemical abundances of 
stars in the Milky Way and M3 1 will be used to infer the past assembly histories of both spiral galaxies as well 
as the structure of their dark matter halos. Complementing the goals of the Gaia mission (V < 17), radial 
velocities and metallicities will be secured for 10 6 Galactic stars to 17< V < 20. Data for fainter stars to 
V ~21 will be secured in areas containing Galactic tidal streams. The M31 campaign will target red giant 
branch stars with 21 < V <22.5 over an unprecedented area of 65 deg 2 . For the extragalactic program, our 
simulations suggest the wide wavelength range of PFS will be particularly powerful in probing the galaxy 
population and its clustering over a wide redshift range and we propose to conduct a color-selected survey of 
1 < z < 2 galaxies and AGN over 16 deg 2 to J ~23.4, yielding a fair sample of galaxies with stellar masses 
above ~ 10 Mg at z ~ 2. A two-tiered survey of higher redshift Lyman break galaxies and Lyman alpha 
emitters will quantify the properties of early systems close to the reionization epoch. PFS will also provide 
unique spectroscopic opportunities beyond these currently-envisaged surveys, particularly in the era of Euclid, 
LSST and TMT 

Subject headings: PFS — cosmology — galactic archaeology — galaxy evolution 



1. INTRODUCTION 

There is currently a major expansion in survey imaging ca- 
pability via the use of CCD and near-infrared detector mo- 
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saics on a wide range of ground-based telescopes. Such imag- 
ing surveys provide accurate photometric and other data to 
enable the study of gravitational lensing signals which trace 
the distribution of dark matter and to conduct census studies 
of Galactic structures and distant star-forming galaxies. For 
over a decade it has been recognized that a similar revolution 
would be provided by a massively-multiplexed spectrograph 
on a large aperture telescope. Spectra provide precise radial 
velocities, metallicities and emission line properties for faint 
and distant sources and enable additional probes of cosmol- 
ogy. The main challenge in realizing this second revolution 
has been access to a wide field telescope, essential for effi- 
cient multi-object spectroscopy of panoramic fields, and the 
cost of implementing the appropriate instrumentation. 

The proposal to construct a Subaru Prime Focus Spectro- 
graph (PFS) emerged following the cancellation in May 2009 
of the Gemini-sponsored Wide-Field Multi-Object Spectrom- 
eter (WFMOS). WFMOS was envisaged as a facility instru- 
ment on the Subaru telescope sharing the optics designed for 
the new prime focus camera, Hyper Suprime-Cam (HSC). 
Two teams received Gemini funding for a conceptual design 
study of WFMOS and, prior to cancellation, a team led by 
Caltech and the Jet Propulsion Laboratory (JPL) secured pre- 
liminary approval. Soon after, however, the Gemini Board 
indicated they did not have sufficient funding to proceed and 
the WFMOS project was terminated. 

The Kavli Institute for the Physics and Mathematics of the 
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Universe (Kavli IPMU) at the University of Tokyo submitted 
a proposal for stimulus funding to the Japanese government 
in September 2009 using design concepts pioneered in the 
WFMOS study led by Caltech and JPL. The successful out- 
come of this proposal in early 2010 initiated the present PFS 
partnership which now includes Caltech/JPL, Princeton and 
Johns Hopkins Universities, the Laboratoire d' Astrophysique 
de Marseille, Academia Sinica Institute of Astronomy & As- 
trophysics (ASIAA) Taiwan, the University of Sao Paulo and 
the Laboratorio Nacional de Astroffsica in Brazil. Our collab- 
oration structure is summarized in Appendix A. 

In addition to the leadership funding provided by Kavli 
IPMU, three important milestones have enabled progress and 
led to the decision to commence construction. Firstly, in Jan- 
uary 201 1 the Subaru Users Meeting endorsed the PFS project 
as a next-generation instrument for the Subaru Prime Focus, 
recognizing the international PFS team and its responsibili- 
ties. This decision led to the establishment of a PFS project 
office at Kavli IPMU in early 201 1 and an allocation of funds 
and manpower by the National Astronomical Observatory of 
Japan (NAOJ) towards integration, commissioning and survey 
operations. A second milestone followed the MOU in Decem- 
ber 201 1 between the Director-General of NAOJ and the Di- 
rector of Kavli IPMU, that anticipates a Subaru Strategic Pro- 
gram for PFS providing up to 300 nights of observing time for 
the PFS team in collaboration with the Japanese astronomical 
community. These developments provided the essential impe- 
tus for the science plans and technical requirements defined in 
this article. The third milestone was a successful Conceptual 
Design Review (CoDR) held in March 2012 which triggered 
the decision to commence construction with first light pre- 
dicted in 2017. The CoDR documentation included a detailed 
science case for PFS and a list of technical requirements and 
this article reproduces these for general interest. 

This article is not intended to provide a technical descrip- 
tion of PFS but a brief overview is helpful (see Fig. 1). Fur- 
ther technical details of the instrument and its current design 
can be found in http://sumire.ipmu.jp/en/2652. PFS is de- 
signed to allow simultaneous low and intermediate-resolution 
spectroscopy of 2400 astronomical targets over a 1.3 degree 
hexagonal field. It shares the Wide Field Corrector and as- 
sociated Field Rotator and Hexapod already constructed for 
the HSC. An array of 2400 optical fibers is in the Prime Fo- 
cus Instrument and each fiber tip position is controlled in- 
plane by a two-stage piezo-electric Fiber Positioner Cobra 
system. Each fiber can be positioned within a particular pa- 
trol region such that these patrol regions fully sample the 1 .3 
deg field. A Fiber Connector relays light to four identical 
fixed-format 3-arm twin-dichroic all-Schmidt Spectrographs 
providing continuous wavelength coverage from 380nm to 
1.3/^m. The blue and red channels will use two Hamamatsu 
2Kx4K edge-buttable fully-depleted CCDs (as in HSC). The 
near-infrared channel will use a new Teledyne 4RG 4K x 4K 
HgCdTe 1.7;um cut-off array. 

The present article describes the detailed scientific case for 
PFS in the context of a Subaru Strategic Program (SSP) of 
~300 nights of observing time. Since such a program would 
not be implemented until 2017 at the earliest, the main moti- 
vation in formulating the team's plans at this stage is in pro- 
viding a list of key science requirements for the technical de- 
sign. In Sections 2 - 4 we describe 3 key components of the 
science case for PFS that will likely form the basis of the Sub- 
aru Strategic Program. For each of these 3 cases we provide 
a summary, a detailed science justification and survey strat- 



egy as well as the flow-down from these to the technical re- 
quirements for the instrument. We summarize these science 
requirements for PFS in more detail in Section 5 and discuss 
some outstanding issues in Section 6 which we will address 
prior to our Preliminary Design Review. 

2. COSMOLOGY 

Summary: PFS will be remarkably powerful in spectro- 
scopic surveys of faint galaxies because of its large multiplex 
gain and the 8.2 meter aperture of the Subaru telescope. The 
extended wavelength coverage provided by the red and near- 
infrared spectrograph arms (650 - 1300nm) will permit a 
unique survey of [O 11] emission-line galaxies extending over 
the redshift range 0.8 < z < 2.4. As large-scale structures 
remain in the linear regime at high redshift, such a survey will 
give detailed new information on the cosmological parame- 
ters as well as the growth rate of structure formation. This 
combination will provide a valuable test of alternative mod- 
els of gravity on large scale which may provide a possible ex- 
planation for dark energy. Multi-color data planned to arrive 
from the HSC imager will be used to select target galaxies for 
spectroscopy and the expected high throughput should yield a 
c^75% success rate of detecting [O lljemission at S/N > 8.5. 
Herein, we propose to conduct a 100 night cosmological sur- 
vey over 1400 deg 2 , sampling galaxies within a comoving vol- 
ume of 9 (Gpc/ft.) 3 over 0.8 < z < 2.4. This will complement 
the lower redshift survey being undertaken by the SDSS BOSS 
collaboration. 

The primary goals of the PFS cosmology survey are to: (1) 
measure the Hubble expansion rate and the angular diame- 
ter distance to 3% fractional accuracies in each of 6 redshift 
bins over 0.8 < z < 2.4 via the baryonic acoustic oscillation 
(BAO) method, (2) use the distance measurements for deter- 
mining the dark energy density parameter ttdc(z) to about 7% 
accuracy in each redshift bin, when combined with lower red- 
shift BAO measurements, (3) use the geometrical constraints 
to determine the curvature parameter to 0.3% accuracy, 
and (4) measure the redshift- space distortion (RSD) in or- 
der to reconstruct the growth rate of large-scale structure to 
6% accuracy since a redshift z = 2.4. These PFS measure- 
ments of the large scale galaxy distribution can be combined 
with complementary weak lensing information from the HSC 
survey in order to significantly improve the cosmological and 
structure growth constraints and reduce uncertainties arising 
from galaxy bias and nonlinear effects that are otherwise ma- 
jor sources of systematic error in spectroscopic surveys. 

2.1. Cosmology Objectives 

The accelerated expansion of the Universe is the most in- 
triguing problem in cosmology. It either requires the introduc- 
tion of a mysterious form of energy, "dark energy", or it could 
signal a breakdown of Einstein theory of General Relativity 
on cosmological scales. To distinguish between these and 
other possibilities requires precise observational constraints 
on both the expansion history of the universe and the growth 
rate of large-scale structure. 

Measurements of galaxy clustering statistics are one of the 
most powerful means of addressing the nature of dark energy. 
The tight coupling between baryons and photons prior to the 
decoupling epoch of z ~ 1100 leaves a characteristic imprint 
on the pattern of galaxy clustering on large scales - the so- 
called baryonic acoustic oscillation (BAO) scale. As the BAO 
length scale is precisely constrained as ~ 150 Mpc from cos- 
mic microwave background (CMB) experiments (Komatsu 
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Fig. 1 . — A brief overview of PFS instruments, which consist of components Wide Field Corrector, Field Rotator, Prime Focus Unit, and Fiber Positioner. 
A Fiber Connector relays light to four identical fixed-format 3-arm twin-dichroic all-Schmidt Spectrographs providing continuous wavelength coverage from 
380nm to 1.3/im. 



et al. 2011), it offers a standard ruler by which we can in- 
fer the angular diameter distance and the Hubble expansion 
rate from the observed correlation function of the galaxy dis- 
tribution. The BAO scale is in the linear or weakly-nonlinear 
density regime and thus provides a robust geometrical test. 
Furthermore, if uncertainties arising from galaxy bias can be 
removed or accurately modeled, we can use the amplitude and 
shape information of the galaxy correlation function in order 
to constrain cosmological parameters as well as the growth 
rate of structure formation. 

Recognizing this, the main scientific questions we seek to 
address with the PFS cosmology survey are: 

1. Is the cosmic acceleration caused by dark energy or 
does it represent a failure of Einstein 's theory of gravity 
on cosmological length scales? 

2. What is the physics of the early universe that generates 
the primordial fluctuations as the seed of large-scale 
structures? 

To address these fundamental questions, the main goals for 
the PFS cosmology survey are to: 

• Constrain the angular diameter distance and the Hubble 
expansion rate via the BAO experiment to the precision 
comparable with, or better than, existing, ongoing or 
planned BAO surveys. 

• Derive the BAO constraints in a redshift range that 
is complementary to those probed by the existing or 



planned BAO surveys on the time scale of the PFS sur- 
vey. 

• Utilize the unique capabilities of the 8.2m Subaai Tele- 
scope and the PFS spectrograph for maximizing cos- 
mological science. 

• Use the shape and amplitude of galaxy correlation func- 
tion in order to constrain cosmological parameters as 
well as the growth rate of structure formation. 

• Combine the weak lensing information, delivered from 
the HSC survey, with the PFS cosmology survey in or- 
der to improve the cosmological constraints by calibrat- 
ing systematic uncertainties that cannot be resolved by 
either of the PFS and HSC surveys alone. 

2.2. PFS Cosmology Survey 

Here we describe the parameters of the PFS cosmology sur- 
vey that are required to meet the above scientific goals. 

Firstly, we will consider which type of galaxy to target with 
PFS. Given the optical and near infrared wavelength coverage 
of PFS, [O II] emission-line galaxies (ELG; [O II] =3121 k) 
are particularly useful tracers allowing an efficient survey out 
to high redshift beyond z = 1, a redshift range that is diffi- 
cult to probe with 4m-class telescopes. Luminous red galaxies 
(LRGs) are a further potentially-useful tracer of large-scale 
structure as studied by the SDSS survey, but at z ^ 1.4 they 
reveal weaker spectral features and are less abundant per unit 
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TABLE 1 
Instrumentation parameters 



Number of fibers 


2400 (600 for each spectrograph) 


Field of view 


1.3 deg (hexagonal - diameter of circumscribed circle) 


Field of view area 


1.098 deg 2 






Fiber diameter 


1.13" diameter at the field center; 1.03" at the edge 




Blue arm 


Red arm 


IR arm 


Wavelength coverage |nmj 


380-670 


650-1000 


970-1300 


Spectral resolution A/ A A 


1900 


2400 


3500 


Pixel scale [A/pix] 


0.71 


0.85 


0.81 


Read-out noise [e — rms/pix] 


3 


3 


4a 


Detector type/read-out mode 


CCD 


CCD 


HgCdTe/SUTR 


Thermal background [e~/pix/sec] 


None 


None 


0.013 


Dark current [e /pix/sec] 


3.89 X 10~ 4 


3.89 x 10~ 4 


0.01 


Spectrograph image quality [fim rms/axis] 


14 


14 b 


14 


Sky continuum 


21.56 mag AB/arcsec 2 @ 1 /xm 


at zenith c 


OH line brightness 


16.6 mag AB/arcsec @ J band at zenith 


Moonlight 


None (dark time) 




Atmospheric extinction 


Variable; continuum is 0.05 mag/airmass @ A > 0.9um c 


Instrument throughput 


Provided from the project office 




Grating wings 


Lorentzian with g of the true number of lines 


Diffuse stray light 


2% of the total light reaching the detector 


Sky subtraction residuals 


2% per pixel d 







NOTE. — PFS instrumentation parameters and associated assumptions used for estimating an expected 
signal-to-noise ratio for an observation of emission-line galaxies. 

11 Per sub-exposure; if the NIR channel is not reset between exposures, the actual assumption corresponds 
to 4\/2 = 5.6 e~ rms. 

b Defocus within the thick CCD is considered separately. 
c Option 0x11003 in the Exposure Time Calculator, v 1 . 
d Equivalent to 1 % per 4-pixel resolution element. 



volume. In the following we focus on ELGs to explore an 
optimal survey design. We may retain LRGs in future consid- 
erations of our survey plans, but their study is not required to 
meet the PFS cosmology objectives defined above. 

2.2.1. Sensitivity of the PFS spectrograph 

To estimate the feasibility of PFS for a wide-field survey of 
ELGs, we have studied the expected performance of measur- 
ing a [O II] line of a galaxy in our targeted redshift range for 
a representative exposure time during the dark lunar phase. In 
doing so, we properly account for the sky emission (contin- 
uum plus OH emission lines) and absorption as well as the 
instrumentation parameters for the current baseline design (as 
listed in Table 1). 

Fig. 2 shows the expected signal-to-noise ratio of the [O II] 
line as a function of redshift, measured with each of the blue, 
red and near-infrared (NIR) arms of PFS. As a working ex- 
ample, here we assume /[on] = 5 x 10~ 17 crg/s/cm 2 for 
the total flux of the [On] doublet, 15min for the exposure 
time, 0.8" seeing size and E(B - V) = 0.05 for the Galac- 
tic dust extinction, respectively. The galaxy radial profile is 
assumed to be an exponential disk with a half-light radius of 
0.3". Note that for the galaxy yield forecasts, we use half- 
light radii from the COSMOS Mock Catalog (Jouvel et al. 
2011), and re-compute the fiber aperture correction for each 
galaxy. 16 We have assumed the 15 minute integration is split 
into 2 sub-exposures for cosmic ray (CR) detection in the 
CCD channel. The NIR channel will perform CR rejection 
by processing of the frames acquired during sample-up-the- 
ramp (SUTR) mode. The cosmology ETC assumes 4 e~ read 
noise per sub-exposure (appropriate for ~ 90 samples along a 

16 This assumes that the [O II] emission traces the /-band continuum in 
which the galaxy sizes were measured. 



450sec ramp). We will probably not reset the NIR channel in 
between sub-exposures, so an overall read noise of AV2 = 5.6 
e~ per pixel would enable the survey described herein. 

In addition to throughput and sky brightness considerations, 
we have considered several other potential limitations. Their 
amplitude is difficult to estimate, but they have been impor- 
tant for previous spectrographs and so we make an explicit al- 
lowance for them so as to adopt a conservative approach. The 
systematic sky subtraction residuals and small-angle stray 
light are very important factors in the study of galaxy spec- 
tra where [O II] is partially blended with a sky line. Diffuse 
stray light is a concern when [O II] lies in a cleaner part of the 
NIR spectrum. 

• Systematic sky subtraction residuals - These are mod- 
eled by adding a "noise" term corresponding to some 
percentage of the sky counts in each spectral pixel. We 
currently set this to 2% of the brightest of the pixel and 
its neighbor on either side (equivalent to 1 % sky sub- 
traction accuracy on a 4-pixel resolution element). 

• Small-angle stray light - We assign to the grating an ef- 
fective number of lines that is 1/3 of the actual number. 

• Diffuse stray light - We take 2% of the OH line flux 
incident on the detector and uniformly spread it over 
all pixels. (This may be appropriate for a detector that 
reflects 10% of the incident radiation, and then there 
are many surfaces that could potentially reflect this ra- 
diation back. Refining this parameter will be a prior- 
ity since the S/N forecasts degrade rapidly if it gets 
worse.) 

Continuing our conservative approach, we assumed the in- 
strumentation throughput at the edge of the focal plane and 26 
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FIG. 2. — Expected signal-to-noise (S/N) ratio for measuring the [O II] emission line as a function of redshift; the blue, green and red curves show the results 
for the PFS blue, red and IR arms in Table 1, respectively, for an total emission line flux of 5 X 10 — 17 erg/cm 2 /s. To properly account for the uncertainties, 
we assumed the instrumentation parameters of the current baseline design listed in Table 1 , an observation at the edge of the focal plane, and included the sky 
emission/absorption and the Galactic dust extinction of E(B — V) = 0.05 and 26 degrees for the zenith angle of the telescope. This computation assumes 15 
min total exposure (split into two exposures; 450sec X 2), cr v = 70 km/s for the velocity dispersion (the intrinsic line width), and 0.8" for the seeing FWHM. We 
also accounted for the finite galaxy size relative to the seeing profile and the fiber size, assuming an exponential profile with scale radius 0.3" for the emission-line 
region. Note that S/N is estimated by the root-sum-square of the spectral pixels (i.e. it is a matched filter combining both doublet members). The current design 
allows a significant detection of [O II] emission line over a wide range of redshift, up to z ~ 2.4 with near-equal sensitivities of the red and NIR arms. 



degrees for the zenith angle. (The latter corresponds to obser- 
vations at declination 5°S, the southern boundary of the HSC 
survey region, and ±0.5 hours away from transit.) 

Fig. 2 shows that the current design of PFS allows a signif- 
icant detection of the [O II] line over a wide range of redshift, 
up to z ~ 2.4. Most importantly, the baseline design provides 
near-equal sensitivity of the red and NIR arms for measuring 
the [O II] line for the same exposure time. Hence PFS can ex- 
ecute a cosmological survey very efficiently over a wide range 
of redshift, provided sufficiently bright ELGs are available for 
study (see below). 

In the following analysis, we set a threshold of S/N = 8.5 
(matched filter) for detection of an ELG. In principle, it may 
be possible to accept less significant detections. However 
given the uncertainties in the airglow and the early stage of 
the instrument design, we consider it prudent to set aside some 
margin in S/N. 

2.2.2. Target selection of emission-line galaxies 

We now address how to optimally select ELGs as suitable 
targets in the proposed redshift range. We will assume that 
we can use the multi-color imaging data of the planned HSC 
survey which will be executed ahead of the PFS survey. The 
currently-planned HSC survey will reach i ~ 26 (5er for a 
point source and 2" aperture), in the 5 passbands grizy over 
~ 1500 square degrees. 

As seen in Fig. 2, if we target ELGs over the wide red- 
shift range 0.8 < z < 2.4, the wide wavelength coverage of 
red and NIR arms allows a very efficient selection of [Oil] 
emission-line galaxies. A g — r color cut is ideal for se- 
lecting galaxies in this redshift range: if an object is blue 



(g — r < 0.3), then it likely has no spectral breaks in the g 
and r bands - this means the redshift is high enough for the 
Balmer/4000A break to have redshifted beyond the r band, 
but the Lyman break has not yet entered the g band. Further- 
more, g — r < 0.3 implies a blue rest-frame UV slope, which 
has a strong correlation with the star-formation activity that 
produces the [O II] emission. 

To estimate the efficiency of various target selections, we 
used the COSMOS Mock Catalog (Jouvel et al. 201 1), where 
fluxes of various emission lines of each galaxy are estimated 
based on physical parameters (SFR, stellar mass and metal- 
licity) using the COSMOS 30 passband photometric data and 
zCOSMOS spectroscopic data. We have chosen the prelimi- 
nary target selection cuts: 

22.8 < g < 24.2 AND - 0.1 < g - r < 0.3 
AND NOT (g > 23.6 AND r — i > 0.3). (1) 

The HSC depth (g, r, i ~ 26 mag AB at 5cr) is sufficient to 
find the target galaxies and to provide accurate g — r and r — 
i colors. The ELGs in the redshift range 0.8 ^ z <, 2.4 are 
primarily selected from the color cut —0.1 < g — r < 0.3, 
the p-magnitude cut gives preference to bright objects while 
reducing low-redshift contamination as can be seen in Fig. 3. 
The condition onr-i for fainter magnitudes is designed to 
tilt the redshift distribution in favor of more objects at z > 1. 
We can further divide the targets into a "bright subsample" 
(g < 23.9) and a "faint subsample" (g > 23.9), with the 
brighter targets prioritized when there is an option to increase 
the success rate. 
In the COSMOS Mock Catalog there are 7847 target galax- 



6 



Objects/deg /cell 



z>0.7 ELG efficiency, f hiz 



1.5 



0.5 



-0.5 




2000 



1000 v 




FIG. 3. — Left panel: The distribution of objects in the COSMOS Mock Catalog on the color-magnitude diagram. Right: The fraction of objects in each cell 
that are z > 0.7 ELGs with [O II] doublets detectable at > 8.5cr in PFS in 2 X 7.5 min exposures. The third dimension (r — i) is not shown on this 2D plot, but 
allows us to select lower or higher redshift galaxies within the PFS survey range. 
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FIG. 4. — The sensitivity of PFS (green curve; 8.5o\ 2 X 450s exposures, 
dark time) to the [O II] doublet at r a g = 0.3" and 1:1 line ratio, versus the 
selected targets (red points). Note that the redshift z = 2.49 corresponds to 
the long wavelength end of the NIR arm. Most targets will yield successful 
redshifts, but some are lost within the atmospheric emission or absorption 
lines, a few are at z > 2.5, and there is a small number of faint blue nearby 
objects (lower-left corner). Note that the line ratio (-F3726 : ^3729) and 
effective radius are re-computed for each galaxy in the COSMOS Mock Cat- 
alog, and hence the sensitivity curve drawn does not correspond to an exact 
boundary between detections and non-detections. 



ies available per PFS field-of-view (1.098 deg 2 for the 1.3° 
FoV diameter). Hence there are a sufficient number of target 
galaxies compared to the number of fibers (iV=2400) for the 
baseline design. The blue-solid curve then shows the redshift 
distribution of [O II] galaxies with S/N > 8.5 for a 15 min- 
utes exposure. To estimate the expected S/N for each galaxy, 
we employed the same method used in Fig. 2 and also used 
the galaxy size information and [O II] doublet ratio available 
from the COSMOS mock catalog. 



Assuming 2400 fibers in the focal plane as in Table 1 and 
using the results in Fig. 5, we can estimate a success rate of 
finding [O II] emission-line galaxies among the target galax- 
ies. We show the magnitude and [O II] flux distributions of 
the targets in Fig. 4, and the redshift histogram of the suc- 
cessful [O II] detections in Fig. 5. Of the targets, 74% have 
successful [O II] detections, 6% are faint blue local (z < 0.3) 
objects, 6% fail due to an [O II] feature that falls off the red 
end of the spectrograph, and the remaining 14% fail due to 
some combination of too faint [O II] feature or overlap with 
an atmospheric emission or absorption complex. To have a 
sufficiently dense sampling of galaxies to trace large-scale 
structures in each redshift slice, we will need multiple vis- 
its of each field; our BAO forecasts found that 2 visits gave 
the best constraints. 

To obtain a reliable estimate of the number of observable 
targets, we took into account the fiber allocation efficiency 
assuming a Poisson distribution of target galaxies on the sky, 
which should be a good approximation for a given wide red- 
shift coverage. We conservatively assume non-overlapping 
patrol zones between the different fibers. 17 The fiber assign- 
ment algorithm is designed to put the easier targets in one of 
the visits ("Visit A") and then the harder targets that require 
better conditions in another ("Visit B"). We divide our targets 
into two tiers - the bright (g < 23.9) and faint (<? > 23.9) sub- 
samples. The fiber assignment logic within each patrol zone 
is then: 

• If at least two bright targets are available, one is as- 
signed to Visit A and another to Visit B. 

• If one bright target and at least one faint target are avail- 

17 That is, we allow a galaxy only to be observed by the nearest fiber, even 
though 21% of galaxies are in overlaps and could potentially be observed by 
either of two fibers. 
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FIG. 5. — The distribution of successful redshifts ([On] detected at 
S/N > 8.5) for the proposed PFS cosmology survey, including breakdown 
into the two visits. The jagged features in the curves reflect the effect of sam- 
pling variance of large-scale structures in the COSMOS field due to the finite 
survey area (the mock is based on the data of 1.24 square degrees). We refer 
to the two visits as "Visit A" and "Visit B" (see text for details), respectively, 
where we preferentially select brighter targets with g < 23.9 in Visit A in 
order to have some flexibility between dark/grey nights. 
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FIG. 6. — An illustration of the BAO features in the power spectrum. The 
cross symbols with error bars show the N-body simulation results, where the 
error bars are the expected la statistical errors of power spectrum measure- 
ment at each fe's bins for the PFS BAO survey of z ~ 1 slices. The solid 
curve shows the analytical prediction at z = I, computed using the refined 
perturbation theory of nonlinear structure formation (Taruya et al. 2009). The 
perturbation theory prediction is in remarkable agreement with the simulation 
result on scales up to k <; 0.25 /iMpc -1 . For comparison, the dotted curve 
represents the prediction from linear theory. 



able, then the bright target is observed in Visit A and 
one of the faint targets is observed in Visit B. 

• If one bright target and no faint targets are available, 
then the bright target is observed in Visit A and the fiber 
becomes a sky fiber in Visit B. 

• If no bright targets are available and there are at least 
two faint targets, then a faint target is observed in both 
Visits A and B. 

• If no bright targets are available and there is only one 
faint target, then the fiber becomes a sky fiber in Visit 
A and the faint target is observed in Visit B. 

• If no targets are available, then the fiber is a sky fiber in 
both Visits A and B. 

This algorithm produces a roughly balanced abundance of sky 
fibers in the two visits. The predicted allocations of the fibers 
are: 

• Visit A: 85% bright targets, 6% faint targets, 9% sky. 

• Visit B: 56% bright targets, 33% faint targets, 11% sky. 

Thus this leave about 240 sky fibers in each visit for a calibra- 
tion of the sky spectrum. 

Note that the two visits could be scheduled in either order. 
Since Visit A has the brighter targets, it can achieve a high 
success fraction under worse conditions than Visit B. We have 
therefore assumed that Visit B takes place during dark time, 
whereas Visit A is scheduled on a night of 7 days from the 
new moon (but at least 45° from the targets). The exposure 
times in both cases are kept at 2 x 450sec. The predicted 
redshift success rate for a S/N > 8.5 threshold is 75% (Visit 
A) 18 or 73% (Visit B). 

18 This would rise to 79% if Visit A were during dark time. 



An alternative to the baseline (using gray time for Visit A) 
would be to use dark time only for the cosmology survey, and 
shorten the exposure time for Visit A, thereby reducing the 
total number of nights required but using time that may be 
in high demand by other programs. This trade will be made 
when we design an integrated observing schedule for PFS (see 
Section 6). 

2.2.3. Survey Strategy 

Using the results of target selection in Fig. 5, we have 
adopted parameters for the PFS cosmology survey summa- 
rized in Table 2. Since our primary observable is the galaxy 
two-point correlation function or the galaxy power spectrum, 
the key factors that govern the results are the geometrical vol- 
ume surveyed and the ratio of clustering power to shot noise, 
UgPg. To have a galaxy power spectrum measurement that 
reaches the sampling variance limit for our volume and is 
not degraded by shot noise, the number density of galaxies 
must satisfy n g P g £ 1 at BAO scales (see Fig. 6). As given 
by the columns of Table 2, the PFS survey we are proposing 
has figPg Z> a few and slightly less than 1 at k = 0.1 and 
0.2 h/Mpc, respectively, over an entire target redshift range. 
With only one visit per field, these numbers are about a factor 
2 smaller than in Table 2. On the other hand, if we have more 
than two visits, the survey area we can cover for a given num- 
ber of nights becomes smaller. Incorporating multiple visits 
ensures more flexibility in optimizing the survey, for exam- 
ple in including a mixture of targets in different magnitude 
ranges. 

We have assumed that the bias factor for the ELGs is given 
by b g (z) = 0.9 + OAz. This specific function was a fit to 
semi-analytic models (Orsi et al. 2010), but compares very 
well to real data: e.g. the DEEP2 "main blue" sample has a 
measured bias of b g = 1.28 ± 0.04 at z = 0.9 (Coil et al. 
2008). 19 We have much less information about clustering of 

19 This was normalized to a erg = 0.9 model; with the lower erg now 
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TABLE 2 

PFS Cosmology Survey Parameters 





Vsurvcy 


JV 9 


n g 


bias 


n g p g( k ) 


n gPg{ k ) 


redshift 


[Gpc/h] 3 


per field 


[10- 4 (h/Mpc) 3 ] 


b g k 


= 0.1/i/Mpc 


k = 0.2/i/Mpc 


0.6 < z < 0.8 


0.59 


85 


1.9 


1.18 


0.74 


0.25 


0.8 < z < 1.0 


0.79 


358 


6.0 


1.26 


2.23 


0.74 


1.0 < z < 1.2 


0.96 


420 


5.S 


1.34 


2.10 


0.68 


1.2 < z < 1.4 


1.09 


640 


7.8 


1.42 


2.64 


0.87 


1.4 < 2 < 1.6 


1.19 


491 


5.5 


1.50 


1.78 


0.59 


1.6 < z < 2.0 


2.58 


598 


3.1 


1.62 


0.95 


0.31 


2.0 < z < 2.4 


2.71 


539 


2.7 


1.78 


0.76 


0.25 



NOTE. — The leftmost column shows the redshift range of each slice, and the other columns show the 
comoving volume (V), the number of [Oil] galaxies per field (Ng), the mean comoving number density 
(fig), the linear bias parameter (b g ) and the values of n g Pg(k) at k = 0.1 and 0.2/i/Mpc for each slice, 
respectively. The survey volume is for a survey area of 1464 square degrees, which is estimated assuming 
15 min of open-shutter time per visit, 2 visits per field, 3 min overhead per visit, and 100 clear nights. 
For comparison, the BOSS BAO survey has the survey parameters: 10000 sq. degrees area coverage over 
0.4 < z < 0.7, Kmrvoy = 4.4 (Gpc//i) 3 , fig = 3 x 10- 4 (/i/Mpc) 3 , b g = 2.3 and n g P g {k = 
0.1 h/Mpc) ~ 5. 



ELGs at redshifts beyond the DEEP2 survey, but Ha emitters 
at 2 = 2.23 have a correlation length of>o = 3.1±0.7Mpc//i 
(Sobral et al. 2010), implying a bias of b g = 1.87t°it 

The power spectrum measurement accuracy depends also 
on the area coverage. In order for the PFS survey to have a 
constraining power of cosmological parameters comparable 
with the existing or planned BAO surveys, we need a suffi- 
ciently large area coverage. We have found that, if about 100 
clear nights are allocated to the PFS cosmology survey, it can 
meet our scientific goals. Hence we assume 100 clear nights 
for the following analysis, and the total area covered is esti- 
mated as 

100 [nights] x 8 [hours] x 60 [min] 
2 [visits] x (15 [min] + 3 [min]) 

x 1.098 [sq. degrees per FoV] = 1464 sq. degreesj(2) 

Here we conservatively assumed 3min overhead for each 
new pointing, which covers readout, slewing and the time 
for accurate fiber positioning, and assumed that 8 hours per 
night are available for observation on source. The comov- 
ing volume in each redshift slice is given in Table 2. The 
total volume is about 9.3 (Gpc//i) 3 , a factor 2 larger vol- 
ume than the SDSS BOSS galaxy survey, which is about 
4.4 (Gpc/h) 3 . A notable strength of the PFS survey is that 
it probes large-scale structure in higher redshifts, where the 
fluctuations are largely in the linear regime and therefore 
allow a cleaner estimation of cosmological parameters. In 
fact the genuine cosmological power comes from the vol- 
ume in Fourier space; the effective volume at each wavenum- 
beris given as V cfi (k) = [n g P g (k)/(l + n g P g {k))f V survcy , 
where V^rvey is the comoving volume. The total number of 
the Fourier modes usable of constraining cosmology is esti- 
mated by integrating the effective volume in Fourier space 
up to the maximum wavenumber k max which is determined 
such that the theoretical model to be compared with the mea- 
surement is reliable up to the /c max . The larger the redshift, 
the higher the wavenumber we can use, because the nonlinear 
scale becomes smaller (higher k). Hence, the proposed PFS 
survey offers much more than a factor 2 improvement com- 
pared to the BOSS constraint (e.g., Anderson et al. 2012). 
Although Table 2 also gives estimates for the lower redshift 
slice 0.6 < z < 0.8, which partially overlaps the ongoing 

favored, the bias would be higher, i.e. "better" for the BAO analysis. 



BOSS and WiggleZ surveys, the cosmological constraining 
power of this slice is not as great due to the smaller areal 
coverage and reduced galaxy number density. However, we 
consider it important to retain this redshift slice as it can give 
a useful benchmark in comparison with other surveys such as 
the BOSS and WiggleZ surveys, particularly for calibrating 
systematic issues. 

2.2.4. Expected cosmological constraints 

We now can estimate the power of the PFS cosmology sur- 
vey in Table 2 for constraining cosmological parameters. To 
ensure a fair comparison with other surveys, we primarily as- 
sess the power of PFS survey in terms of its BAO geometrical 
constraints. 

Geometrical constraints: 

The galaxy two-point correlations are measured as a function 
of the separation lengths between paired galaxies. The po- 
sition of each galaxy needs to be inferred from the measured 
redshift and angular position. Then the separation lengths per- 
pendicular and parallel to the line-of-sight direction from the 
measured quantities are given as r± oc Ad and rii oc Az, 
where A9 and Az are the differences between the angular po- 
sitions and the redshifts of the paired galaxies. For this con- 
version, we need to assume a reference cosmological model 
to relate the observables (AO, Az) to the quantities (r±,ru). 
Thus, the wavenumbers are given as 



k± 



ref 



Da{z) 

DA,rei( z ) 



k±, 



, grefjf) , 

Kll.rcf = „/ * fell 



H(z) 



(3) 



The quantities with subscript "ref" are the quantities esti- 
mated from the observables assuming a "reference" cosmo- 
logical model, and the quantities without the subscript are 
the underlying true values. Since the reference cosmolog- 
ical model assumed generally differs from the underlying 
true cosmology, it causes an apparent distortion in the two- 
dimensional pattern of galaxy clustering. In principle, the 
distortion could be measured using only the isotropy of clus- 
tering statistics (Alcock & Paczynski 1979), but a more ro- 
bust measurement can be obtained using features in the power 
spectrum, particularly if they are at a known scale so that we 
can measure both Da(z) and H(z). In particular, the CMB- 
inferred BAO scale of 150 Mpc gives a powerful standard 
ruler for this geometrical test (Eisenstein et al. 2005; Perci- 
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val et al. 2007; Blake et al. 201 1; Hu & Haiman 2003; Seo & 
Eisenstein 2003). 

The use of this scale without edge effects in a survey field 
requires the survey to be contiguous on a scale large compared 
to the BAO length; at our minimum redshift (z = 0.8), 2.5 
BAO lengths corresponds to 7.5 degrees on the sky, so we set 
this as our minimum width. This requirement will be refined 
further by simulations. 

In more detail, the galaxy power spectrum in redshift space 
is given in the linear regime as 



-Pff,s(fc_L,rcf,fc||,rof;2 : ) 



D A ^ cl {z) 2 H{z) 
H Iet {z)D A (z) 2 



xb 2 P^(k;z) + P sn 



(4) 



where b g is the linear bias parameter, j3 is the linear 
redshift-space distortion (RSD) parameter, defined as (3 = 
(l/bg) dhxD/d\na\ z (Kaiser 1987), D is the linear growth 
rate, P„ is the linear mass power spectrum, and P sn is a pa- 
rameter to model the residual shot noise. We can use the BAO 
features in the linear power spectrum P^ n as a standard ruler 
in order to constrain Da(z) and H(z). The BAO constraints 
are relatively robust against the galaxy bias uncertainty and 
the other nonlinearity effects, because none of the systematic 
effects introduces any particular length scale comparable with 
the BAO scale. Further, if we can use the shape and amplitude 
information in the galaxy power spectrum, we can constrain 
the growth rate as well as other cosmological parameters such 
as the neutrino mass and the primordial power spectrum pa- 
rameters (Takada et al. 2006), as we will discuss below. 

To make the parameter forecast, we have used the method 
developed in Seo & Eisenstein (2007). In this method, we in- 
clude the smearing effect of the BAO features due to the bulk 
flow of galaxies in large-scale structure (Matsubara 2008; 
Taruya et al. 2009; Nishimichi & Taruya 201 1). For the BAO 
survey of multiple redshift bins, the Fisher information matrix 
of model parameters can be computed as 



1 r k 



- 2Trk 2 dk din P SjS (fc, [i\ z±) 



x , femln 2(2^)3 d Pa 
din P g , s {k, fi; Zi) 



X ■ 



-V cf t(k;zi 



x exp -fc^Si - fcV(£| - Si) , (5) 

where /i is the cosine between the wavevector and the line- 
of-sight direction, /i = fcy/fc; ^ z is the sum over different 

redshift bins; dP g s /dp a is the partial derivative of the galaxy 
power spectrum (Eq. 4) with respect to the a-th parameter 
around the fiducial cosmological model; the effective survey 
volume V c ff and the Lagrangian displacement fields Em and £ 
to model the smearing effect are given as 



V e g(k, n;zi) 



n g (zi)Pg, s (k,ij;Zi) + 1 
Ei(z)sc rcc D(z)So, 

E||(«) = C r ecD(z)(l + / fl )E . 



(*)(6) 

(7) 
(8) 



Here V survcy (zi) is the comoving volume of the redshift slice 
centered at zf, the present-day Lagrangian displacement field 
is E = ll/i -1 Mpc for er 8 = 0.8 (Eisenstein et al. 2007); 



D(z) is the growth rate normalized as D(z = 0) = 1; 
f g = din D '/din a. The parameter c roc is a parameter to 
model the reconstruction method of the BAO peaks (see be- 
low). In Eq. (5), we take the exponential factor of the smear- 
ing effect outside of the derivatives of P giS . This is equiva- 
lent to marginalizing over uncertainties in £|| and Ej_. The 
growth rate in E || or E^ takes into account the smaller smear- 
ing effect at higher redshift due to the less evolved large-scale 
structure. For the parameters, we included the cosmological 
parameters, the distances in each redshift slice, and the nui- 
sance parameters: 

Pa = {^mo, A s ,n s , a s ,tt m0 h 2 , fl h oh 2 , D A (zi), H(z i ) 1 

6 5 (z,),/3(^),^n(^)}, (9) 

where A s , n s and a s are parameters of the primordial power 
spectrum; A s is the amplitude of the primordial curvature 
perturbation, and n s and a s are the spectral tilt and the run- 
ning spectral index. The set of cosmological parameters de- 
termines the shape of the linear power spectrum. By using 
the method above, we can estimate the cosmological distance 
information solely from the BAO peaks, not from the broad- 
band shape of the power spectrum. For the fc-integration, we 
set fc m in = 10~ 4 /i/Mpc and fc max = 0.5 h/Mpc for all the 
redshift slices, but the exponential factor in Eq. (5) suppresses 
the information from the nonlinear scales. The Fisher param- 
eter forecasts depend on the fiducial cosmological model for 
which we assumed the model consistent with the WMAP 7- 
year data (Komatsu et al. 201 1). 

Further, we assume that we can implement the promising 
method of Eisenstein et al. (2007) for improving the BAO 
measurements. Since the peculiar velocity field of galaxies in 
large-scale structure can be inferred from the measured galaxy 
distribution, the inferred velocity field allows for pulling back 
each galaxy to its position at an earlier epoch and then recon- 
structing the galaxy distribution more in the linear regime. As 
a result, one can correct to some extent the smearing effect 
in Eq. (5) and sharpen the BAO peaks in the galaxy power 
spectrum. Recently, Padmanabhan et al. (2012) implemented 
this method with real data from the SDSS DR7 LRG catalog, 
and showed that the reconstruction method can improve the 
distance error by a factor of 2. The improvement was equiv- 
alent to reducing the nonlinear smoothing scale from 8.1 to 
E n i = 4.4 h~ 1 Mpc, about a factor 2 reduction in the displace- 
ment field. To implement this reconstruction method requires 
a sufficiently high number density of the sampled galaxies in 
order to reliably infer the peculiar velocity field from the mea- 
sured galaxy distribution. Each redshift slice of the PFS sur- 
vey (see Table 2) satisfies the requirement; the number density 
of galaxies in each redshift slice is higher than that of both the 
SDSS DR7 LRGs (n ~ 10" 4 (/i/Mpc) 3 ) and the BOSS LRGs 
(h ~ 3 x 10~ 4 (/i/Mpc) 3 ). Hence we can safely assume that 
the reconstruction method can be applied to the PFS BAO ex- 
periment. In the Fisher matrix calculation, we used c roc = 0.5 
for an implementation of the reconstruction method 20 . 

Finally, we have used the CMB information expected from 
the Planck satellite experiment, which gives a precise con- 
straint on the sound horizon scale in order for us to use the 
BAO scale as a standard ruler. In addition to the cosmo- 
logical parameters (Sl m o , A s , n s ,a s , f2 m o/i 2 , tlboh 2 ), we in- 
cluded, in the CMB Fisher matrix, r (the optical depth to the 

20 Then we found we can roughly reproduce the distance measurement 
accuracy for the SDSS LRGs as found in Padmanabhan et al. (2012). 
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last scattering surface) and Da(zcmb) (the angular diameter 
distance to the last scattering surface). Then we can com- 
pute the Fisher matrix for the BAO experiment by adding the 
galaxy and CMB Fisher matrices; a , , = F^ xy + F^ ck . 
The dimension of the Fisher matrix for the PFS survey in 
combination with the expected Planck information is 38 x 38 
(see Eq. 9). When further combined with the SDSS and 
BOSS BAO information, the dimension of the Fisher matrix 
increases accordingly. 

Fig. 7 shows the expected accuracies of determining the an- 
gular diameter distance and the Hubble expansion rate in each 
redshift slice with the PFS cosmology survey (Table 2). The 
errors include marginalization over uncertainties of other pa- 
rameters. The PFS forecasts can be compared with the accu- 
racies of the existing and ongoing SDSS/BOSS surveys. As 
can be clearly seen, the PFS cosmology survey can constrain 
Da(z) and H(z) over a wider range of redshift yet with sim- 
ilar precision to the SDSS and BOSS surveys. Even though 
the PFS area coverage is smaller than that of SDSS or BOSS 
surveys by a factor of 10 (1400 vs. 10000 sq. degrees), the 
PFS survey covers a factor 10 and 2 times larger comoving 
volume than the SDSS and BOSS surveys, respectively. 

Cosmological implications: 

It is worth noting that the BOSS survey will give the cluster- 
ing statistics of the Lyman-a forest over 2.1 ;$ z ^ 3.4. PFS 
thus serves to fill a natural 'gap' in-between the galaxy and 
Lyman-a BAO experiments, allowing for probing the expan- 



sion history over the entire range of redshifts, 



3, i.e. 



through the period where it is believed that the cosmic ex- 
pansion went from decelerated to accelerated phases. If we 
model the expansion history as parametrized by the dark en- 
ergy model (wq, w a ) (Linder 2003) and the curvature param- 
eter (K), 



H\z) = H* 



+ ^dc,0 a 



-3(l+w +w a ) e 3w a (a-l) ^0) 



we can propagate the distance measurement errors into the 
accuracies of estimating the parameters. To be more explicit, 
we can do this, based on the Fisher information matrix for- 
malism, by projecting the BAO Fisher matrix into different 
parameter space: 



F n 'h 



d Pa psub d Pb 
-fab 



dp a > 



dpb 



(ID 



Here the new set of parameters p a > is given as p a > = 
(fideo, fimo'' 2 ! , Wq, w a ), which specifies the cosmic ex- 
pansion history given by the above equation, and F^ h 
is the sub-matrix computed by inverting the sub-matrix 
of the inverse of the full BAO matrix, [F _1 ] Q( g, contain- 
ing only the parts of the geometrical parameters, p a = 
{ilmo, fi m o/i 2 , D^Zi), H(zi)}. Hence the derived con- 
straints on p a ' include marginalization over other parameters 
such as the galaxy bias and the (5 parameters. Table 3 shows 
the expected accuracies of the dark energy parameters and the 
curvature parameter for the PFS survey. Here w p i vo t is the 
dark energy equation state at the "pivot" redshift, at which the 
dark energy equation of state is best constrained for a given 
survey. The quantity FoMdc is the dark energy figure-of-merit 
defined in the Dark Energy Task Force Report (Albrecht et al. 



2006), which quantifies the ability of a given survey for con- 
straining both w and w a ; FoM dc = l/[<T(wp ivot )(j(w a )}, 
which is proportional to the area of the marginalized con- 
straint ellipse in a sub-space of (wo,w a )- Table 3 clearly 
shows that the PFS BAO can significantly tighten the param- 
eter constraints over the SDSS and BOSS surveys. Most in- 
terestingly, the PFS has the potential to constrain the curva- 
ture parameter to a precision of 0.3%. If we can detect a 
non-zero curvature, this would represent a fundamental dis- 
covery giving a critical constraint on the physics of the early 
universe, for example insight into different inflation scenarios 
(Efstathiou 2003; Contaldi et al. 2003; Freivogel et al. 2006; 
Kleban & Schillo 2012; Guth & Nomura 2012). 

Nature of dark energy: 

The parametrization (wo,w a ) adopted for the dark energy 
equation of state samples only a narrow range of dark en- 
ergy models. Given there is no well-accepted model for dark 
energy we seek to interpret our PFS data in a more model- 
independent way. The wide redshift coverage of PFS surveys, 
in combination with the SDSS and BOSS surveys, allows us 
to do this by directly reconstructing the dark energy density 
as a function of redshift solely based on the geometrical BAO 
constraints. To illustrate this, we use the Hubble expansion 
history parametrized in terms of dark energy density parame- 
ters in each redshift bins: 



Hi 



Ono(l + Z) 3 



Hi 



(l+^) 2 + 



Pde,z,( z G z i) 



PcrO 



(12) 

where pd c ,zi is the dark energy parameter in the redshift bin 
centered at Zj. For the combined BAO survey of SDSS, BOSS 
and PFS, we include the 9 dark energy densities, pde{zi), 
given in 9 redshift bins (8 redshift bins of the galaxy sur- 
veys plus the redshift bin from 2.4 to Zcmb)- Then, similarly 
to the method described around Eq. (11), we can propagate 
the BAO-measured distance errors into the accuracies of re- 
constructing the dark energy densities in each redshift bin. 
Fig. 8 shows the result, where we assumed the cosmological 
constant, pd c = constant, as the fiducial model. The figure 
clearly shows that the PFS BAO survey is capable of recon- 
structing dark energy densities up to high redshift, z ~ 2 for 
a model in which pdc(z) — constant, thereby testing various 
types of early dark energy models. 

Further, we should emphasize that adding the PFS BAO 
measurements to the SDSS and BOSS information can 
significantly improve the reconstruction of dark energy 
density at low redshift, partly because the SDSS+BOSS 
BAO alone cannot break degeneracies between parameters 
(f2 m o, f! m o, K, pde(zi)) given the distance measurements to 
the three redshifts (zsdsSj z bosSj zcmb), an d a lso because 
the angular diameter distances for the PFS redshifts are all 
sensitive to dark energy densities at low redshift via the inte- 
gration relation between the angular diameter distance and the 
Hubble expansion rate. However, we should note that the con- 
straint on the curvature parameter is significantly degraded in 
this case to cf(VIk) = 0.057 from the result in Table 3 because 
of the larger number of free parameters. Thus the accuracy 
of the curvature parameter is sensitive to which dark energy 
model we use for marginalization. 

Testing the growth of structure: 

We now turn to the utility of measuring the redshift- space dis- 
tortion (RSD) effect and the broad-band shape of the galaxy 
power spectrum. If we can reliably model the RSD and the 
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FIG. 7. — Fractional errors of determining the angular diameter distance and the Hubble expansion rate via the PFS BAO experiment (see Table 2) including 
marginalization over uncertainties of other parameters. The expected accuracies are compared to the existing and ongoing SDSS and BOSS surveys. The PFS 
survey will provide geometrical constraints to higher redshift than the SDSS and BOSS surveys, but with comparable precision. The solid curve in each panel 
shows the fractional difference when changing the dark energy equation of state parameter from the fiducial model to = — 1 to w = —0.9. 



TABLE 3 

Forecasted accuracies of cosmolooical parameters 



Survey 


S 2dc0 


Wpivot 


W a 


FoMde 


(Ik 


m„,tot [eV] 


/nl 


SDSS+BOSS 


0.0061 


0.076 


1.2 


11 


0.0071 


0.188 


16 


SDSS+BOSS+PFS 


0.0051 


0.059 


0.36 


47 


0.0030 


0.133 


11 



NOTE. — The constraints on fidco, w p i vot , w a and Qk(= —K/Hq) are from the 
BAO distance measurements in Fig. 7, i.e. not including the information on the broad- 
band shape of the galaxy power spectrum. Note that ui p i vo t is the dark energy equation 
state at the "pivot" redshift, at which the dark energy equation state parameter is best 
constrained for the given PFS BAO measurements. The constraints on the neutrino 
mass m„,tot and /nl are derived by including the broad-band shape information. See 
text for the details. 



shape of the power spectrum in the weakly nonlinear regime 
including a possible scale-dependent bias, we can use this in- 
formation not only to improve the cosmological constraints 
(Takada et al. 2006), but also to constrain the growth rate 
which is sensitive to theory of gravity on cosmological scales. 
Encouraging progress is being made via many efforts to de- 
velop a more accurate model of the redshift-space power spec- 
trum in the weakly nonlinear regime (see Fig. 6) (Matsubara 
2008; Taruya et al. 2009; Nishimichi & Taruya 2011; Tang 
et al. 2011; Hikage et al. 2012). 

To estimate the power of the PFS survey, we use the lin- 
ear theory prediction for the amplitude of the RSD effect, 
(3(z) = f g (z)/bg(z), in Eq. (4), where f g is defined by the 
growth rate as f g = din D/ In a. Then we can include the 
RSD effect in the Fisher matrix formalism by using f g in each 
redshift slice instead of treating f3 as parameters (see Eqs. 4 



- 9). With this implementation, we can break degeneracies 
between the RSD effect f g /b g and the galaxy bias uncer- 
tainty b g from the measured anisotropic modulations in the 
redshift-space galaxy power spectrum. Then we can in turn 
use the amplitude and shape information of the underlying 
linear power spectrum. 

Fig. 9 shows the expected accuracies of constraining the 
growth rate, f g {= dlnD/d\na), in each redshift slice via 
the RSD measurements. The PFS survey can constrain the 
growth rate in each redshift to a 6% accuracy. In particu- 
lar, PFS will provide accurate constraints on the growth rate 
at redshifts beyond z = 1, when the cosmic expansion is in 
its decelerated phase. Such constraint are very important for 
testing whether dark energy is an illusion caused by an incom- 
plete understanding of General Relativity. 
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Reconstruction of dark energy density at each redshfit 
SDSS+BOSS 



SDSS+BOSS+PFS 
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redshift: z 

FIG. 8. — Expected accuracy of reconstructing the dark energy density pa- 
rameter at each redshift, Qdc(z) = Pdc(z)/[3.ff 2 0z)/87rG] from the BAO- 
measured Da(z) and H(z) in Fig. 7. Here we considered the cosmological 
constant (pdc( z ) = PdoO =constant) and the fiat universe (fi^ = 0) as 
the fiducial model. Adding the PFS BAO constraints to the SDSS and BOSS 
constraints enables reconstruction of the dark energy density to z ~ 2, and 
also significantly improves the precision at low redshifts, as the comoving 
distance at the high redshift arises from an integration of H(z). 
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FIG. 9. — Marginalized errors of reconstructing the growth rate, f g 
d\nD/ In a, in each redshift slice. 



Other constraints: 

With the growth rate constraints and the information on the 
shape of the galaxy power spectrum, we can also constrain 
other interesting parameters such as the sum of neutrino 
masses (m^tot) and the degree of primordial non-Gaussianity 
(/nl). Primordial non-Gaussianity induces a characteristic 
scale-dependent biasing effect on the galaxy distribution at 
very large scales (Dalai et al. 2008) that are well in the linear 
regime and cannot be explained by other nonlinearity effects. 
Hence we can use the largest-scale signal of galaxy cluster- 
ing to explore the signature of the primordial non-Gaussianity. 
Table 3 shows the expected accuracy of constraining /nl to an 
accuracy of <t(/nl) — H- PFS does not have a comparable 
power of constraining /nl to the constraints expected from 
other probes such as the Planck experiment (/nl ~ 5) due to 
the relatively small area coverage, which limits the access to 
the largest-length scales. 



On the other hands, the neutrinos of finite mass scale, 
as found by terrestrial experiments, cause a suppression 
in the galaxy clustering power on scales smaller than the 
neutrino free-streaming scale, which imprints a character- 
istic scale-dependent effect on the galaxy power spectrum 
(Takada et al. 2006). The amount of the suppression ef- 
fect scales with the sum of neutrino mass as AP g /P g ~ 
— 80,/o/Omo — — 8m I/i tot/(94.1 eVfi m n/i 2 ); the neutrinos of 
m^tot =0.1 eV, close to the lower bound of the inverted 
neutrino mass hierarchy, leads to about 6% suppression in 
the galaxy power spectrum compared to the case without the 
massive neutrinos. Hence, we can use the measured cluster- 
ing amplitude to constrain the neutrino mass. However, the 
achievable precision of neutrino mass depends on the level of 
our understanding on the nonlinear power spectrum including 
the galaxy bias uncertainty (Saito et al. 2008, 2009). Here, by 
assuming that an accurate model of the galaxy power spec- 
trum is available, we estimate the power of PFS for constrain- 
ing the neutrino mass. To be more precise, we assumed that 
the following set of parameters, instead of Eq. (9), can model 
the redshift-space galaxy power spectrum based on the ex- 
tended perturbation theory based method in combination with 
numerical simulations: 

p a = { O m0 , A s , n s , a s , Q m0 h 2 , fl h0 h 2 , fl K , w , w a , 

m v , tot, b g {zi),P S n{zi)}. (13) 

In this parameter estimation we did not use the reconstruction 
method (i.e., set c rec = 1 for the reconstruction parameter in 
Eq. 8), because the reconstruction method of BAO peaks al- 
ters the shape and amplitude of power spectrum. With this 
implementation, we can include the shape and amplitude in- 
formation of the power spectrum for constraining the cosmo- 
logical parameters, marginalized over uncertainties of the nui- 
sance parameters. As can be found from Table 3, the PFS 
survey can achieve a precision of <7(m„ t ot) —0.13 eV. 

Synergy with HSC: 

Finally we note there will be many synergistic opportunities 
enabled by the fact that the PFS survey will be undertaken in 
the same sky areas as the HSC imaging survey. Weak lens- 
ing lensing information from HSC will be very effective in 
correcting and calibrating systematic effects inherent in the 
galaxy clustering analysis, the nonlinear redshift-space dis- 
tortion and the galaxy bias uncertainty, up to the z ~ 2 slice 
(Hikage et al. 2012; Nishizawa et al. 2012). The spectro- 
scopic data from PFS survey can likewise be used to calibrate 
the photo-z errors or the redshift distribution of HSC imag- 
ing galaxies, which is one of the major uncertainties in the 
HSC cosmology. Thus, by combining the HSC imaging and 
PFS spectroscopic surveys, we can significantly improve the 
cosmological constraints making the joint HSC and PFS ex- 
periments comparable to a Stage-IV Dark Energy experiment 
in the parlance of recent US studies. 

2.3. Scientific Requirements for PFS Cosmology Survey 

As discussed above, PFS has the unique capability to ex- 
ecute a very powerful cosmology survey across a wide range 
of redshifts, extending considerably current and planned BAO 
surveys on < 4 m-class telescopes. Here we summarize the 
survey and instrumental requirements for the cosmological 
applications. These requirements are listed in Tables 4 and 
5. We regard these as minimum requirements and, of course, 
where possible will strive to improve upon these subject to 
cost and schedule impact. 
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TABLE 4 

Cosmology Survey Requirements 



Science yield requirements 


Distance measurements 

Dark energy reconstruction 

Curvature 

Growth of structure 


^ 3% measurement of Da (z) and H (z) in each of 6 redshift bins 
via BAO (0.8-1.0, 1.0-1.2, 1.2-1.4, 1.4-1.6, 1.6-2.0, and 2.0-2.4) 
;$ 7% measurement of Q^ c (z) in each of 6 redshift bins via BAO 
Measure U K to < 0.3% via BAO 

;$ 6% measurement of the growth rate of structure in each of 6 bins via RSD 


Galaxy catalog requirements 


Redshift range 

Number density of galaxies 

dN/dz of ELGs 

Total survey area 

Incorrect redshift fraction 

Redshift precision, accuracy 

Survey geometry 


0.8 < z < 2.4 (0.8 < z < 1.6 minimum) 

> 2900 deg~ 2 

n g p 9 > 1 (0-8 < z < 1.6) or figP g > 0.5 (1.6 < z < 2.4) at k = 0.1/i/Mpc 

> 1400 deg 2 
< 1% 

Az/(1 + z) < 0.0007, la (a v < 200 km/s) 

Width > 7.5 degrees; < 4 contiguously-connected survey regions 


Survey implementation requirements 


Total nights 
Lunar phase 
Imaging survey 


~ 100 clear nights 

Dark ( 1 of 2 visits) or age < 7 days (other visit) 
HSC gri data to « 26th magnitude AB (5<r) 





TABLE 5 




PFS Instrument Requirements 


Wavelength coverage 


650 < A < 1300nm (650 < A < lOOOnm minimum) 


Number of fibers 


2400 


Overhead 


< 0.2 X the open-shutter time 


Throughput 


Average > 22% (red) or > 24% (NIR) 




Worst part of band, > 20% (red) or > 18% (NIR) 




(Excludes atmosphere, central obscuration + WFC vignetting, and fiber aperture 




effect.) 


Fiber aperture factor 


Encircled energy in fiber is > 59% (point source) or > 45% (galaxy, r e g = 0.3") 




(Equivalent to 0.8" FWHM seeing + 11 /im rms/axis additional aberration 




+ 0.12" fiber offset.) 


Spectrograph image quality 


< 14 fim rms per axis (excluding fiber geometric size, pixel tophat, and 




internal defocus due to thickness of red CCD) 


Spectral resolution 


Red: R ~ 3000; NIR: R ~ 4000 (to resolve out OH lines but limit read noise) 


Stray light 


Near OH lines: Lorentzian wings at < 3x amplitude of perfect grating 




Diffuse: equivalent to < 2% of total sky brightness spread over detector 


Read noise 


< 3 (red) or < 4 (NIR) e~ rms per pixel 




(If the NIR channel is not reset between exposures, < 4V2 e~ rms is acceptable.) 


Sky subtraction accuracy 


< 1% of sky background per 4-pixel resolution element 



3. GALACTIC ARCHAEOLOGY: NEAR-FIELD 
COSMOLOGY 

Summary: The nature of dark matter and its role in galaxy 
formation are best tested on the scale of or below, an indi- 
vidual galaxy. We propose a PFS survey of ancient resolved 
stars in the Milky Way and in the Andromeda Galaxy (M31) 
to address this question. There is significant tension between 
the predictions of the standard Cold Dark Matter cosmogony 
and the observations of galaxies on small scales, such that 
one must appeal either to a baryonic or dark matter solution. 
As the baryonic physics governing how galaxies fanned and 
evolved can be traced by the properties of stars of all ages, 
including those from the earliest epochs, PFS can contribute 
fundamentally to this area of near-field cosmology. 

Our proposed PFS Galactic Archaeology ( GA) Survey aims 
to measure radial velocities and metallicities for a large sam- 
ple of stars, in synergy with both HSC and the ESA space mis- 
sion Gaia. It can be viewed as an extension of SDSS/SEGUE, 
going significantly deeper. The targets consist of ancient 
Galactic stars belonging to the stellar halo and disks, both 
in the smooth components and in streams and satellites, and 
bright red giants in the halo of M31. Our proposed survey 
makes full use of the current spectrograph design, i.e. re- 
solving power of R = 2, 000 to 3,000, wavelength cover- 



age of 3,800 A to 9,000 A with 2,000-3,000 fibers per PFS 
field. Importantly, it has already been demonstrated by the 
SDSS/SEGUE team that stellar spectra we will obtain with 
this resolution can provide the necessary high quality data 
for our science. 

3.1. Galactic Archeology Objectives 

Our understanding of how galaxies like the Milky Way and 
M3 1 formed and evolved in the expanding Universe remains 
enigmatic. The current paradigm of cosmic structure forma- 
tion based on Cold Dark Matter (CDM) models proposes that 
galaxy formation is driven by hierarchical assembly of dark 
matter halos, starting from sub-galactic scales. The repeated 
merging and overall clustering of many small dark halos over 
cosmic time forms a larger, galaxy-sized halo. The associated 
cooling and collapse of baryonic matter confined in the dark 
halos at each stage of the hierarchy forms the visible galaxies, 
with the morphology of the final galaxy set by the physics of 
star formation and of the merging process. Indeed, mergers 
are predicted to play a key role in the creation of stellar halos, 
bulges, thick disks and even some parts of thin disks. 

However, CDM models have encountered several fun- 
damental difficulties in explaining observations on galactic 
scales consistent with an excess of small scale power. One 
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FIG. 10. — An illustration of how the number of objects of a given type de- 
pends on their present day mass for different dark matter models (Ostriker & 
Steinhardt 2003) including Collisionless Cold dark matter (CCDM), Strongly 
Self-Interacting dark matter (SIDM), Warm dark matter (WDM), Repulsive 
dark matter (RDM), Fuzzy dark matter (FDM), Self-Annihilating dark mat- 
ter (SADM), Decaying dark matter (DDM) and Massive Black Holes (BH). 
Observations on small scales are key to distinguishing these. 

of the most serious issues is the prediction of many more sub- 
halos in a Milky Way-sized system than the modest number 
of visible satellites. This discrepancy can only be reconciled 
if the vast majority of small halos are dark. They would then 
only reveal their presence through their dynamical effects on 
visible systems. Our understanding of the star formation pro- 
cess on these small scales, which correspond to dwarf satel- 
lite galaxies, remains limited, so there may well be a large 
population of dark satellites. Alternatively, the assumption 
of collisionless CDM models may not be correct (Ostriker & 
Steinhardt 2003, see Fig. 10 for the predictions of various dark 
matter models). The merging history of the 'building blocks' 
of galaxies is determined by the power spectrum, i.e. by the 
nature of dark matter, and this we plan to investigate through 
the study of the stellar populations in both the Milky Way and 
M31. 

Our primary science goal with Subaru/PFS is thus to test the 
assembly process of CDM on Galactic scales, by dedicated 
observations of old stars nearby, formed at high look-back 
times. Resolved ancient stars in the Milky Way and other 
Local Group galaxies are ideal targets because they offer us 
our most detailed views of galactic structure and evolution 
through their kinematics and chemical abundances (Freeman 
& Bland-Hawthorn 2002). Kinematics of stars reflect past 
galaxy collapse and/or merging events and their distribution 
in phase space as defined by integrals of motion such as an- 
gular momentum remains basically unchanged (Fig. 11). The 
chemical abundances of stars reflect their past star-formation 
history and chemical evolution, possibly in association with 
the dynamical state of proto-galactic clouds such as their col- 
lision and merging. All of these processes are actually con- 
trolled by hierarchical clustering of CDM via self-gravity on 
galactic and sub-galactic scales. Thus, to assess what CDM 
models predict, it is essential to derive the spatial distribu- 
tion of dark matter subhalos in the Galaxy and Andromeda 
through their dynamical effects on visible stellar systems. 
Stars are indeed ideal tracers of a background gravitational 
field dominated by dark matter. 
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FIG. 11. — Left: Model distribution of tidal streams in spatial coordi- 
nates where the different colors represent different satellites (from Freeman 
& Bland-Hawthorn 2002). These stream-like features disappear after several 
dynamical times. Right: Model distribution of nearby stars in the integrals of 
motion space, i.e., E vs. L z and E vs. L, based on numerical simulations 
of falling satellites into the Milky Way (from Helmi & Zeeuw 2000). The 
different colors represent different satellites. Shown is the final distribution 
of stars after 12 Gyr within about 6 kpc from the Sun, after convolution with 
the errors expected for Gaia. It is clear that each of the progenitor galaxies 
can be traced via the current phase space distribution. 



To make progress we propose to obtain spectra for about 
a million stars in the Milky Way and Andromeda, as one 
component within the framework of the proposed PFS Sub- 
aru Strategic Program. These spectroscopic studies will be in 
perfect synergy with the upcoming astrometric satellite, Gaia, 
that will provide very precise measurements of distances and 
proper motions for Galactic stars, and the HSC survey which 
will provide the essential large imaging data set. Therefore, 
PFS, in combination with Gaia and HSC, will provide us with 
unique data in the area of near-field cosmology. With these 
data we hope to gain an ultimate understanding of the nature 
of dark matter and the associated galaxy formation processes. 

The main questions we seek to address in our dedicated 
Galactic Archaeology (GA) survey are summarized as 
follows. 



1. What is the merging history of the Milky Way? - ad- 
dressing the role and nature of dark matter in galaxy 
formation. 

2. How did the old Galactic components (thick disk and 
stellar halo) form? - addressing baryonic physics in 
the context of CDM. 
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3. How does M31 differ from the Milky Way? - contrast 
merging and baryonic processes on small scales in two 
systems. 

3.2. PFS Galactic Archaeology Survey 

A summary of our GA survey plan is shown in Table 6. 
Its primary goal is to measure radial velocities and overall 
metallicities for about a million stars belonging to the field 
stellar halo and thick disk in the Milky Way, together with 
stars associated with substructure - streams and satellites - 
and bright red giant stars in a wide area of M31's halo, along 
its minor axis. Our proposed survey will allow metallicities to 
be determined to — 0.2 dex using the SDSS/SEGUE stellar- 
parameter pipeline (as detailed below) and radial velocities to 
a few km s _1 . The survey comprises two parts. 

3.2. 1 . The Milky Way Survey 

The primary targets in this survey lie in the old stellar com- 
ponents of the Milky Way. PFS will measure spectra of stars 
with 17 < V < 21.5 mag, thereby complementing data from 
Gaia. The Gaia mission will provide astrometric observa- 
tions for ~ 1 billion stars with a precision of 12-25 p,as for 
V < 15 and 300 p,as for V < 20, whereas radial velocities 
and metal abundances will only be measured for bright stars 
with V < 17 and V < 12, respectively. The Milky Way sur- 
vey is therefore primarily aimed at obtaining radial velocities 
and metal abundances for two samples. In the first Gaia sam- 
ple, we will select stars in the magnitude range of V < 20 
mag, for many of which accurate distances and full spatial 
motions will be available by the time of our PFS survey and 
the full final data will be released by 2022. Targeting F/G 
dwarf stars is especially important to densely map the phase 
space distribution and find kinematic streams (Fig. 1 1) within 
distances of ~ 10 kpc from the Sun. In the second Faint sam- 
ple, PFS will target stars beyond the Gaia astrometry limit, 
i.e., with V > 20 and as faint as V = 21.5 to map the dense 
phase-space distribution up to 30 kpc and to reach very dis- 
tant, giant stars at more than 100 kpc. This sample is selected 
from stars associated with halo substructures in the form of 
stellar streams (Sagittarius Stream and Orphan Stream) and 
those in the outskirts of the closest ultra-faint Galactic dSph 
(Segue I and associated 300 km s" 1 stream). This sample 
also targets general field stars located in the outer-most parts 
of the disk and halo from beyond the range of Gaia (V > 20) 
to the PFS limit, which corresponds to a distance range from 
60 to 130 kpc for a luminous giant, and from 100 to 200 kpc 
for a HB or RR Lyr star. Measurements of radial velocities as 
well as metal abundances of these stars will set important con- 
straints on chemo-dynamical properties of the outer reaches of 
the disk and halo. 

This PFS survey will complement the on-going Gaia-ESO 
Survey with VLT in the Southern hemisphere. Gaia-ESO is 
a dedicated spectroscopic survey, led by GGilmore, target- 
ing > 10 5 stars, systematically covering the various struc- 
tural components of the Galaxy, providing homogeneous data 
of radial velocities and elemental abundances also in synergy 
with Gaia astrometry. Based on high-resolution spectra taken 
from two wavelength windows (R = 19, 800 for 5339-5619 
A and R = 16, 200 for 8484-9000 A) and 110 science fibers, 
the survey is measuring multi-element chemical abundances 
of stars over 300 nights. The team has made significant ef- 
forts to define spectroscopic standard fields, some of which 
are equatorial and so can be accessed from Mauna Kea. 



We now discuss the two components of our PFS Milky Way 
survey in more detail: 

( i) Gaia sample: 

The justification for studying around a million stars in syn- 
ergy with Gaia is as follows. For mapping of the phase-space 
distribution for different stellar populations, we assume that 
the stellar masses of the building blocks for the halo/thick disk 
range from 10 3 to 10 7 M Q as deduced from the luminosities 
of Galactic satellites (e.g., Belokurov et al. 2006). Then the 
stellar halo whose total mass is ~ 10 9 M Q is comprised of 
about ~ 10 4 building blocks (assuming a power law mass 
function with index in the range —1.8 to —1.5). If at least 100 
member stars are needed for each building block to charac- 
terize its velocity and metallicity distribution with 10 % sta- 
tistical significance, we need to observe ~ 1 million stars in 
total. 

For separating and characterizing each of the Galactic com- 
ponents through its velocity, space and metallicity distribu- 
tions, we require at least 40 velocity bins for —200 to 200 
km s _1 with Av = 10 km s , 100 spatial positions with 10 
bins for each of two dimensional (R, z), (b, I), or (along, per- 
pendicular to) the line of sight, and 2 to 3 metallicity ranges. 
Again, provided that at least 100 stars at each velocity bin are 
needed to characterize its velocity distribution with 10 % sta- 
tistical significance, we require a sample of a million stars in 
total. 

This science is most robustly achieved with F/G-type dwarf 
stars at distances up to 10 kpc, i.e., consisting of both halo 
and disk stars in roughly equal proportions. Based on ac- 
curate space motions and overall metallicities available for 
these stars, phase-space and metallicity distributions of each 
of the stellar systems/substructures can be derived. As de- 
tailed below (see Fig. 16), we will have about 300 F/G dwarfs 
at I = 180° per PFS field per 2 magnitude interval in the rel- 
evant magnitude range, suggesting about 500 F/G dwarfs per 
PFS field for 17 < V < 20. To enable ancillary science using 
other stellar populations, no color cut will be adopted for the 
Gaia sample. Many of these additional subsamples will con- 
tain stars in the boundary of inner/outer halo at around r ~ 20 
kpc whose full spatial motions will be determined for the first 
time. Such data will provide an excellent tracer of the under- 
lying gravitational potential dominated by the dark halo and 
hence the density distribution of dark matter in its inner part 
of r < 20 kpc (e.g., Holmberg & Flynn 2000; Garbari et al. 
2011; Smith et al. 2012). 

To acquire the necessary data for as many as 1 million stars, 
the total survey area must be at least ~ 400 sq. degree (~ 300 
PFS pointings) at b ~ 30 deg, in which ~ 55 %, ~ 40 % 
and ~ 5 % belong to the halo, thick disk and thin disk, re- 
spectively. We require spectra with a S/N~ 30 per resolution 
element at A = 5000 A for precisions in metal abundance to 
0.2 dex and radial velocity to a few km s _1 , on the basis of 
the SDSS/SEGUE calibration method as detailed below. 

( ii) Faint sample: 

The Milky Way halo contains numerous tidally elongated 
stellar streams covering the sky, as depicted in Fig. 12 (Be- 
lokurov et al. 2006). The most impressive feature is the stream 
emanating from the Sagittarius dSph galaxy together with a 
second, fainter stream nearby. The Sgr dSph galaxy has been 
hailed as the 'Rosetta Stone' for the merging processes since 
its discovery in 1994 (Ibata et al. 1994) but we still do not 
understand its detailed origin (for example, how many orbits 
it has been losing material) or the kinematics and metallic- 
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TABLE 6 

Summary of Galactic Archaeology Survey 



Survey 


Magnitude range 
(mag) 


Area 
(deg 2 ) 


No. Fields 


Survey Time 
(nights) 


Comments 


LR-mode survey: 












1. The Milky Way 


17 < V < 21.5 




(208) 




(1) b = 30, < I < 270 








(46) 




(2) b = -30, 60 < I < 120 








(10) 




(3) b = 60, 1 = 90 








(12) 




(4) Outer disk I = 180 








(24) 




(5) 'Field of Streams' 






390 


300 


75 


Sum of (l)-(5) 


2. M31 halo 


21.5 < V < 22.5 


65 


50 


30 


HSC sample 


Total 




455 


350 


105 






FIG. 12. — Stellar substructures revealed in the SDSS imaging data (taken from Belokurov et al. 2006), and proposed PFS pointings (filled red circles) over 
these streams and the surroundings of the faint satellite, Segue I. 



ity distribution of the associated streams. Although the lower 
stream appears to represent tidal material from the Sgr dSph, 
the upper stream cannot be explained currently (e.g., Fell- 
hauer et al. 2006; Penarrubia et al. 2010) (see also Koposov 
et al. 201 1, for the extension of these streams in the Southern 
Galactic hemisphere). The Orphan Stream is another striking 
feature. Although numerous faint dSphs have been revealed 
by SDSS in its general vicinity, its origin remains unclear. 

To better understand the origin and dynamical importance 
of these streams and their association with satellite galaxies, 
or more generally, with building blocks and their assembly 
process of the stellar halo, PFS will undertake a compre- 
hensive chemo-dynamical mapping of the region containing 
the Sgr Stream, Orphan Stream and the outskirts of Segue I 
(Fig. 12). This satellite galaxy shows an extended stellar 
distribution possibly due to tidal interaction and/or overlap- 
ping fore/background stellar populations (see Belokurov et al. 
2008; Simon et al. 2010, for more details). Recently Carlin et 
al. (2012) suggested that Segue I may be associated with the 
tidal disruption of the progenitor galaxy of the Virgo Over- 
density. 

In addition to such faint, known substructures beyond the 
reach of Gaia, large parts of the Galactic components espe- 
cially at the faint limits are still unexplored. Examples in- 
clude the outer parts of the thin disk at I = 180°, where 
Monoceros stream (which is a Northern hemisphere target) 
is located, as yet unidentified kinematically cold substruc- 



tures at high Galactic latitudes at b ~ 60°, the interface of 
halo/disk/bulge components and more. Subaru/PFS will al- 
low us to measure detailed velocity distributions of these faint 
parts of the Galactic components. By the time of the PFS 
survey, proper motions available from ground-based imaging 
survey, such as PanSTARRS, combined with the photometric 
distance obtained from the stellar parameters will provide us 
with the six dimensional phase space distribution at distances 
beyond those reached by Gaia. 

We propose to observe a total area of ~ 390 sq. degree 
(~ 300 pointings) covering all the survey fields for both the 
Gaia and faint samples. Survey fields include stripes along 
constant Galactic latitudes at around b = 30° with 0° < I < 
270° (208 pointings) and b = -30° with 60° < / < 120° 
(46 pointings), high latitude fields (b = 60°) at I = 90° (10 
pointings), the outer-disk fields at I = 180° (12 pointings) and 
the Field of Streams (24 pointings). A total of ~ 75 nights is 
required assuming 2 hr exposures to obtain S/N^ 30 spectra 
for stars with V — 21.5 mag and that 4 fields per night are 
observed. For the brighter stars with V < 20, grey time can 
be used. 

3.2.2. The M31 halo surx'ey 

The primary targets in the M31 survey are bright red gi- 
ants with 21.5 < V < 22.5 mag, i.e., stars around the tip 
of the RGB with / ~ 20.5 selected along the minor axis of 
the galaxy (Fig. 13). Pre-imaging observations and selec- 
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FIG. 13. — Proposed PFS pointings along the minor axis of M31's halo (~ 50 pointings), in which several stellar streams are included. The map of RGB stars 
is taken from Richardson et al. (201 1). 



tion of the candidate halo giants will be provided by an HSC 
imaging survey and we will utilize a newly-developed narrow- 
band filter (NB515 filter with CW = 5145 A and BW = 80 A, 
which is similar to the DD051 filter 21 ), to assist in removing 
foreground Galactic dwarfs. A survey along the minor axis 
will enable us to derive the chemo-dynamical properties of the 
general M31 halo as well as to study in detail several stellar 
streams, including a cold stream-like feature in the north-west 
part of the halo, part of which corresponds to Stream F discov- 
ered by previous Subaru observations (Tanaka et al. 2010). 
The distribution of stars along these M3 1 streams will be es- 
pecially useful in constraining the number of orbiting dark 
matter subhalos, which induce visible density variations along 
the stream through dynamical disturbances (Carlberg 2011). 
A careful subtraction of foreground Milky Way stars through 
spectroscopic observations will be crucial for constraining the 
number of orbiting CDM subhalos in Andromeda and com- 
paring with theoretical predictions. 

PFS will measure radial velocities and metal abundances 
for M31 halo stars to set important constraints on their in- 
trinsic chemo-dynamical properties. Individually-determined 
spectroscopic metallicities will represent a major advance 
over those derived photometrically or by stacking spectra. To 
estimate the number of RGB stars in the M3 1 halo per PFS 
field, we follow the recent observations with Keck/DEIMOS 
by Gilbert et al. (2012), in which the pre-selection of can- 
didate RGB targets using the DD051 filter was made. At 
a projected distance of 80 kpc from the M31 center, we 
expect ~ 2600 objects per PFS field (including significant 
foreground/background contamination even after NB5 15 pre- 

21 The DD051 filter is an intermediate filter in the Washington photomet- 
ric system, where the strong surface gravity sensitivity of the Mg I triplet and 
MgH band features near 5150 A allows discrimination between dwarfs and 
giants of F-K spectral type. 



selection), among which about 100 secure RGB stars of the 
M31 halo can be extracted. We propose to observe ~ 66 sq. 
degrees (~ 50 pointings) as shown in Fig. 13. In order to ob- 
tain spectra with S/N~ 20 per resolution element at A = 5000 
A, 5 hr exposures are required at V = 22.5 mag. Thus, ~30 
nights are required for this component of the survey. 

3.3. Deriving stellar atmospheric parameters from PFS 
spectra 

We now briefly discuss how we will derive the key mea- 
sures of chemical composition essential to our science objec- 
tives from our PFS spectra. 

A number of techniques estimating the stellar metallicity 
[Fe/H] and other atmospheric parameters (effective tempera- 
ture T c g and surface gravity log g) are available to exploit our 
PFS spectra, in combination with photometric data provided 
by, e.g., SDSS and HSC. Well established techniques include, 
e.g., fitting observed spectra with grids of synthetic ones for 
obtaining the most likely set of ([Fe/H], T c g , log g) (e.g., Wil- 
helm et al. 1999; Allende Prieto et al. 2008; Kirby et al. 2008, 
2009, 2010; Re Fiorentin et al. 2007; Kordopatis et al. 201 1), 
using the Ca II K line, auto-correlation analysis of the stellar 
spectrum, and Ca II triplet lines for obtaining [Fe/H] (Beers 
et al. 1999; Norris et al. 2008; Starkenburg et al. 2010), and 
using gravity-sensitive Ca I line and Mg features for obtain- 
ing log g (Morrison et al. 2003). Each of these methods is 
usually optimal over restricted ranges of color (g — r), metal- 
licity [Fe/H], and signal-to-noise ratio S/N. A great advantage 
of PFS in this respect is that PFS stellar spectra are designed 
to cover a broad wavelength range from the blue (~ 3900 A) 
to red (~ 9, 000 A) continuously, so that it will be possible to 
apply various calibration methods and arrive at a robust set of 
stellar parameters ([Fe/H], T c g, logg). 

Using such multiple techniques for the estimation of stel- 
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FIG. 14. — Spectra of metal-poor dwarf stars in the Milky Way with R. = 2, 000, which are obtained by degrading Subara/HDS spectra with R. = 50, 000 
taken from Ishigaki et al. (2012). This spectral region contains dominant a-element spectral features in the wavelength window, which can be used to determine 
[a/Fe] with A[a/Fe] as small as 0.1 dex for S/N~ 30 spectra (Lee et al. 2011). The black and red lines show dwarf stars, BD+13 2995 and G 18—24 with 
[Mg/Fe]= 0.3 and 0.1, respectively. 



lar parameters, the SEGUE Stellar-Parameter Pipeline (SSPP) 
has been developed and utilized in the SDSS/SEGUE sur- 
vey (Lee et al. 2008; Allende Prieto et al. 2008). Recent 
work (e.g., Li et al. 2010) has successfully demonstrated that 
the SSPP can provide reasonable estimation of stellar atmo- 
spheric parameters, even for non-SEGUE data with similar 
spectral resolutions (as low as R = 1,000 ~ 2,000), pro- 
vided the wavelength coverage extends from 3800 A (to in- 
clude Ca II HK) to 5200 A (to include Mg features at 5150 
A which are sensitive to surface gravity). Adopting the SSPP 
(and implements its subsequent improvements up to the time 
of our PFS survey) as the calibration pipeline of PFS spec- 
tra will enable us to obtain, for our sample of a million stars, 
metallicities [Fe/H] to a precision of 0.2 dex. Furthermore, 
with this pipeline, the [a/Fe] ratio commonly used to distin- 
guish enrichment by Type II or la SNe, will be secured to a 
precision of ~ 0.1 to 0.2 dex for spectra with a S/N= 30 (Lee 
et al. 2011). Fig. 14 shows spectra of metal-poor dwarf stars 
in the Milky Way with R = 2, 000, where the spectral region 
contains dominant a-element spectral features, which can be 
used for the calibration of the [a/Fe] ratio, as demonstrated 
by Lee et al. (2011). 

Flux calibrated spectra will provide independent measure- 
ments of stellar atmospheric parameters, and can be advanta- 
geous for the abundance measurements. Regular observations 
of flux standards should be sufficient to define the relative cal- 



ibration to better than 5 %. Absolute flux calibration is not 
necessary since the absolute fluxes can be deriving by adjust- 
ing the relative fluxes to match the extant photometry data. 

3.4. Technical Requirements 

We now summarize the technical requirements essential for 
the GA PFS survey to meet its scientific objectives. 

As discussed above, we require spectra with S/N > 30 per 
resolution element for a broad range of wavelengths, includ- 
ing the Ca II lines (3933 A), in order to determine metallicities 
better than ~ 0.2 dex and radial velocities better than ~ 10 
km s ~ 1 . The latter velocity accuracy is essential to distinguish 
stellar components in velocity or phase space and corresponds 
to a precision better than the velocity dispersion of the Galac- 
tic disk or typical dSphs. The wavelength region around 5150 
A is also important as it includes Mg lb and MgH features 
which provide good sensitivity to surface gravity. 

We need to observe many target stars per PFS field in or- 
der to assemble a statistically meaningful sample. Based on 
a standard star-count model by Robin et al. (2003), there are 
~ 2, 300 stars of V = 17 mag at Galactic latitude b = 30 deg 
per PFS field and a comparable number of stars is expected 
at fainter magnitudes (Fig. 15). We have further examined 
the SDSS DR8 to estimate the observed number of stars per 
PFS FoV at various Galactic longitudes and latitudes within 
SDSS footprints. Fig. 16 illustrates the number of stars with a 
given (extinction corrected) g magnitude range at I = 180° as 
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FIG. 15. — Star count predictions (per magnitude and per deg 2 ) from the 
Besancon model of the Milky Way (Robin et al. 2003) in the V band at I = 
90° or 270°, for latitudes 10° to 90° from top to bottom (0°, 20°, 45° and 
90° with solid lines, 10°, 30° and 60° with dashed lines). For each latitude 
the bifurcation is due to the warp, the highest part being for £ = 90°, b > 
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FIG. 16. — Number of stars from the SDSS DR8 with a given go magnitude 
range per PFS field of view at I = 180° is plotted against Galactic latitudes. 
Left and right panels are for the magnitude ranges of 17.5 < go < 19.5 and 
19.5 < go < 21.5, respectively. 



a function of Galactic latitudes, b. The total number of stars 
per PFS FoV with g < 21.5 is as high as 3000 at b = 30° and 
1500 at b = 60°. We will set no color cuts in selecting targets 
but will utilize all stars for different science aspects. Thus, 
a multiplex gain of 2000-3000 will be sufficient for the pro- 



posed Galactic studies. Also with this setup, it is possible to 
target a sufficient number of stars in the direction of the M3 1 
halo (red giants in M3 1 and foreground/background contami- 
nation) as already mentioned in the previous subsection. 

4. GALAXY EVOLUTION 

Summary: The goal of the PFS Galaxy Evolution Survey is 
to follow the growth of the full panoply of the galaxy popula- 
tion from cosmic dawn to the present. While many of the basic 
physical processes that drive galaxy evolution (dark matter 
halo merging, gas accretion, star formation and associated 
energy release, galaxy merging, black hole accretion and as- 
sociated energy release) are known, how and when they op- 
erate remain unknown. We will use the unprecedented wave- 
length coverage of PFS, in particular its 1 — 1.3/im near- 
infrared coverage, to explore the redshift range 1.4 < z < 2.2 
when the star formation rate density and black hole growth 
were at their peak. Using Ly a emission, we will trace the 
growth of galaxies and black holes all the way to the epoch of 
reionization, 2 < z < 7. With the deep broad- and narrow- 
band imaging from the HSC survey, we will study these young 
galaxies with unprecedented statistics using PFS. We propose 
a 100 night survey covering 16 deg 2 with 3 main compo- 
nents: ( i) a color-selected survey of half a million galaxies 
with 1 < z < 2 to Jab = 23.4 mag, with a z < 1 component 
limited to Jab = 21, (ii) a survey of 30, 000 bright dropout 
galaxies and Ly a emitters over 2 < z < 7 and ( Hi) a survey 
of color-selected quasars from 3 < z < 7. 

4. 1 . Galaxy Evolution Objectives 

SDSS has provided a detailed understanding of the optical 
properties and large-scale distribution of present-day galaxies. 
Fainter surveys such as zCOSMOS (Lilly et al. 2009), DEEP2 
(Cooper et al. 2006), and VVDS (Le Fevre et al. 2005) have 
studied areas of < 2 square degrees at z ~ 1; at this epoch 
galaxies had similar morphologies, luminosity ranges, and en- 
vironmental characteristics as they do today. However, the 
rates per unit comoving volume of star formation and black 
hole growth were at their highest at somewhat higher red- 
shifts 1 < z < 2 (Fig. 17). We seek to understand how 
galaxies, gas, and dark matter halos are distributed at these 
epochs, how each of these components interact (via merging, 
cooling flows, star formation, energy feedback, and reioniza- 
tion), and whether star formation was fundamentally different 
from what we deduce today. To date there are spectroscopic 
samples of only a few thousand blue, star-forming galaxies 
with z > 1.4, covering < 1 deg 2 (e.g., Steidel et al. 2010). 
Particularly difficult to study is the so-called redshift desert 
1.4 < z < 2.2, since no strong and reliable redshift indica- 
tors fall in the optical bandpass. The PFS spectrograph, with a 
wavelength coverage of 0.38-1.3 /im, is specifically designed 
to explore the redshift desert over volumes comparable to the 
SDSS at low redshift. 

PFS can also study galaxies from z ~ 2 to the redshift fron- 
tier of z ~ 7 (Fig. 17). Here we will utilize HSC to target Ly- 
man break galaxies (LBGs) and Ly a emitters (LAEs), obtain- 
ing unprecedented numbers of spectra of both populations. 
PFS will thereby probe the epoch when the star formation of 
early galaxies changed the ionization state of the intergalactic 
medium (IGM). Its panoramic capability will enable studies 
of the topology of ionized bubbles as well as the clustering 
properties of young galaxies at 2 < z < 7. 

It is now recognized that all massive galaxies contain super- 
massive black holes (SMBHs) in their centers, with masses 
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Fig. 17. — Left: Star formation rate per comoving volume as a function of redshift (Hopkins & Beacom 2006). Different symbols show 
measurements made at different wavelengths. The solid lines show parametric fits to the data (refer to Hopkins & Beacom for details). The PFS 
redshift ranges are shown as black dashed vertical lines. Right: Expected number of clusters with M > 5 x 10 13 Mq per deg 2 survey volume 
with the stated depths. The estimates use the "lightcone" simulations of Sehgal et al. (2010). and populate halos with galaxies using the halo 
occupation distribution model from Tinker & Wetzel (2010). 



ranging from 10 5 — 10 10 Mq. The evolution of black hole 
growth with cosmic time resembles that of the star-formation 
activity. The steep drop-off of both at z < 2 points to a link 
between supermassive black hole growth and galaxy forma- 
tion. Moreover, the masses of SMBHs are correlated with 
those of their host bulges in the present-day Universe, sug- 
gesting that they evolved together (e.g., Kormendy & Rich- 
stone 1995; Magorrian et al. 1998; Marconi & Hunt 2003). A 
full understanding of galaxy evolution requires study of the 
accretion history of SMBHs. PFS will target color-selected 
quasars ~ three magnitudes fainter than the SDSS to a redshift 
of z = 7 and will study their clustering properties, the evolu- 
tion of BH mass density with cosmic time, and the evolution 
in the metal content and ionization state of the inter-galactic 
medium. 

The wide wavelength coverage and large area of the PFS 
will enable an enormously powerful galaxy redshift survey. 
We envision a survey of 16 deg 2 over ~ 100 nights. At 

1 < z < 2, we will survey to Jab ~ 23.4 mag with 3 hour 
exposures, yielding a fair sample with stellar masses above 
~ 10 10 Mq at z rj 2. To cover comparable objects at lower 
redshifts, we will also survey z < 1 targets down to J = 21 
mag with 2 hour integrations. We will measure redshifts for 
~ half a million galaxies allowing us to determine luminos- 
ity/mass functions that are not cosmic-variance limited. At 

2 S z ~ 7 we w iH target bright, star-forming galaxies. We 
expect to observe 15,000 LBGs and LAEs, more than an order 
of magnitude more than existing spectroscopic samples. We 
will focus on three imaging surveys for target selection, based 
on the Hyper-Suprime Cam (HSC) imaging survey. The sur- 
vey will have three components; we focus on the HSC Deep 
layer, which will reach a 5 a depth of % < 27.2 mag. We will 
target 16 deg 2 of the HSC Deep survey that contains existing 
deep NIR imaging down to Jab = 23.4 mag. HSC Ultradeep 
covers 3.5 deg 2 and we will perform deeper spectroscopy in 
this area. Finally, HSC Wide covers 1400 deg 2 , over which 
we will search for rare high-redshift quasars (z > 3). 

The scientific objectives of the PFS Galaxy Evolution Sur- 
vey are as follows: 



1. The Build-up of Stellar Mass Density: We will derive 
stellar masses and star formation rates from the spec- 
tra (Fig. 19). It has been shown that stellar population 
modeling works well even at our low resolution (e.g., 
Kriek et al. 2006, 201 1; Panter et al. 2007; Chen et al. 
2012). We will measure star formation rates both from 
the modeling and from [O II] emission out to z ~ 2.5, 
which we will calibrate using Ha at low redshift. At 
yet higher redshift, we can use the UV continuum as 
a star formation rate indicator (Fig. 19). We will trace 
the stellar mass function and star formation history of 
galaxies as a function of color and galaxy density at 
1 < z < 2. 

2. The Growth of Structure: We will measure the spa- 
tial correlation functions of galaxies on both small and 
large scales (i.e., one and two halo terms; Fig. 18 
right). Spectroscopic redshifts give significantly better 
clustering measurements than do photometric redshifts. 
On large scales, we will measure the galaxy bias as a 
function of galaxy properties and redshift, which is a 
key constraint on galaxy formation models. On smaller 
scales, we can study the properties of galaxies in cluster 
and group environments, and see how galaxy formation 
is related to environment and dark matter halo mass. 
We will tie the evolution in stellar mass to that of the 
underlying dark halos from 1 < z < 7. We will also 
study the cluster/proto-cluster population as a function 
of redshift. Over the survey volume, we expect to find 
— 325 clusters with mass M > 5 x 10 13 M (and —65 
massive clusters with M > 10 14 Mq) at 1 < z < 2 
(Fig. 17), an unprecedented sample when compared to 
the very few found in current deep pencil-beam sur- 
veys. At 2 < z < 6 we expect to find — 100 proto- 
clusters using LBGs and LAEs as tracers. 

3. Gas Inflow and Outflow: We will trace the interplay be- 
tween gas accretion and feedback using both direct and 
indirect means. We will trace the mass-metallicity rela- 
tion using strong emission-line diagnostics for z < 1.6 
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Fig. 18. — Left: Depth versus redshift for existing and planned large redshift surveys. The symbol size represents survey area. Each survey 
is placed roughly at the median redshift of the survey. PFS occupies a unique position in this parameter space, as it is the only large survey 
capable of filling in the redshift regime between 1.5 < z < 2. Right: This figure shows the number of spectroscopic pairs separated by less 
than one arcminute. It highlights the power of PFS to (i) study small-scale clustering on the group scale and (ii) probe the gas distribution in 
galaxy halos using absorption line probes. Fiber-based spectroscopic surveys doing only one pass over the sky suffer from the fiber collision 
limitation and cannot access small-scale pairs. We note that we have not counted overlap regions for the BigBOSS survey. 



and with stellar photospheric lines in the UV at z > 2 
by stacking spectra (Ando et al. 2007). We will search 
for outflow or inflow using interstellar absorption lines 
such as Mg n, S IV, C IV, etc in stacked spectra. We 
will also stack background galaxy and quasar spectra 
as a function of impact parameter around foreground 
galaxies to map out the gas kinematics in the outer ha- 
los of star forming galaxies (e.g., Steidel et al. 2010; 
Tumlinson et al. 201 1; Bordoloi et al. 201 1). Using the 
galaxies themselves as a backlight, we will search for 
the signature of cold-gas inflow and test the cold ac- 
cretion hypothesis with a sample ~20 times larger than 
that of Steidel et al. 

4. The Build-up of Supermassive Black Holes: We will 
measure the quasar luminosity function from 3 < z < 
7 to ~ 3 mag deeper than the SDSS. A wide-area 
search of ;>1000 deg 2 is needed to find rare luminous 
z > 6 quasars, while we will study fainter populations 
over the 16 deg 2 galaxy survey. Using quasar emis- 
sion lines, we will study the evolution of the BH mass 
function, and the accretion history of these black holes. 
We will be able to study the clustering properties of 
the quasar population using cross-correlation with the 
dropout galaxies in our sample, which will reveal the 
typical dark halo mass of the quasar population as a 
function of redshift and luminosity (e.g., Adelberger & 
Steidel 2005), and relate the relative growth of galax- 
ies and BHs. Finally, we will uncover obscured and 
low-luminosity active galaxies using emission-line di- 
agnostics in the 1 < z < 2 survey. 

5. The Epoch of Reionization, Ionized bubbles: For the 
first time, we will constrain the topology of ionized 
bubbles at z ~ 7 with a complete spectroscopic sam- 
ple of LAEs. Current surveys (e.g., Ouchi et al. 2010) 
have shown the promise of using the mean Ly a profile 
shape as a tracer of neutrality, but the small samples of- 



fer no information on large-scale spatial variations. An 
ambitious but practical opportunity will be to exploit 
the ~ 9, 000 LAEs from our PFS survey to investigate 
the possibility of detecting evolution in Ly a profiles 
over scales of 10 Mpc, comparable to the expected size 
of typical ionized bubbles. We will calibrate the intrin- 
sic Ly a lineshape using a critical sample of LAEs at 
z ~ 2.3, for which we will have robust systemic red- 
shifts from [Oil]. 

6. The Epoch of Reionization, Neutral Fraction: We will 
utilize the evolution of the Lya luminosity function 
(e.g., Ouchi et al. 2010) and the complementary mea- 
sure of the fraction of galaxies with Ly a emission with 
unprecedented large samples of LAEs and LBGs, re- 
spectively (e.g., Stark et al. 2011). This dual approach 
will allow us calibrate possible intrinsic evolution in the 
Lya luminosity function arising from changes in the co- 
moving star formation density and linking emission line 
demographics to LBG colors will address possible ef- 
fects of dust extinction. The resulting decline in the vis- 
ibility of Ly a over 6 < z < 7 will then be interpreted 
in terms of a change in the neutral fraction using appro- 
priate models (e.g., McQuinn et al. 2007) . We will ob- 
tain complementary constraints on the neutral fraction 
using the Gunn-Peterson troughs of z ~ 6 — 7 quasars. 

The wide wavelength coverage of PFS is an unique asset. 
The NIR arm allows us to measure the [O II] line (and thereby 
reliable redshifts) out to z = 2.5 (Fig. 19). It allows us 
to trace galaxy populations via their diagnostic spectral fea- 
tures from low to high redshifts with no gaps. As illustrated 
in Fig. 19, we will have access to many important emission 
and absorption features over an unprecedented range of red- 
shifts, which will enable new and important science. We will 
cross-calibrate star-formation and metallicity indicators such 
as [O II] and Ha; both can be observed to z ~ 1. 
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FIG. 19. — Redshift range over which various spectroscopic properties of 
galaxies can be tracked assuming a flux-limited survey down to i = 23.5. 
Left: The capabilities of PFS with the full wavelength coverage (opti- 
cal+NIR). Right: The same, but with only the optical band. The dark and 
light areas show the redshift range where a primary and secondary set of 
spectral features will be available for galaxies with masses above a charac- 
teristic value (M*). Passive and star forming galaxies are marked in red and 
blue, respectively. 



4.2. Simulated Data 

We have generated realistic simulated spectra to determine 
what types of measurements are feasible as a function of con- 
tinuum and line luminosity. Our model for the throughput 
of the spectrograph includes the detector quantum efficiency, 
light-losses at the fiber due to seeing, losses in the fiber and 
fiber coupling, reflectivity of the Subaru primary mirror and 
transmission of the corrector, throughput of the camera and 
collimator (including vignetting) and dichroics as a function 
of wavelength, the grating transmission, and the atmospheric 
transparency. It is encouraging that our treatment is similar 
to (but independent of) that described in §2.2.1, and that we 
derive similar results. 

The most important, and most challenging, aspect of our 
simulations involves modeling the sky subtraction. Sky emis- 
sion lines are a major source of noise in the red and NIR arms. 
Again, our treatment is similar to that outlined in §2. We as- 
sume first that we can measure the overall sky level at a given 
time at the 0.5% level. This requirement is theoretically pos- 
sible but technically challenging in practice. Then, for each 
sub-exposure of 20 min, a random-deviate wavelength error 
of 0.01 pixel and a 0.5% amplitude error in the subtracted 
sky, again easily supported by the Poisson noise in the lines 
themselves, is added and propagated forward into the final 
combined spectrum. Finally, we add 1 percent of the total 
flux on the detector as a smooth scattered background. 

The spectra are combined with inverse-variance weight- 
ing to downweight pixels that are contaminated by the sky 
(no clipping is performed), and then resampled at R = 
400, 300, 300 for the blue, red, and NIR arms respectively 
(Figs. 20 and 21) for study of the continuum S/N. When 
we examine emission lines, we instead use a spectral res- 
olution R ~ 1200 to properly resolve emission lines with 



a w 100 kms - . R w 300 is not ideal for measuring spec- 
tral absorption features either. However, because such a large 
fraction of the spectrum is eaten by night sky lines, it is con- 
venient to examine spectra at this resolution as representative 
of the information content in the spectra. 

In this subsection we will show that our main science goals 
are achievable. We assume exposure times of 0.3-3 hrs for 
the continuum-selected survey, spending three hours on the 
sources in the faintest half-magnitude bin (Jab = 23 — 23.4 
mag). Further, we plan to reposition fibers rapidly, thereby 
allowing us to leave bright targets after as short as a single 
20 min integration. We assume six-hour integration times for 
fainter LBGs and LAEs. 

We have two primary sources for the spectral energy dis- 
tributions that we use. For the continuum-selected survey, we 
use model fits to NEWFIRM Medium-Band Survey photome- 
try for galaxies spanning the color-magnitude diagram. These 
fits, from the code EAZY (Brammeret al. 2008) include emis- 
sion lines. In Fig. 20, we show a representative red and blue 
galaxy at our detection limit of Jab = 23.4 mag at z w 2, 
observed in a 3hr integration with PFS. The figure shows that 
the 4000 A break is detected with enough significance to al- 
low redshift detection at the magnitude limit even without the 
presence of emission lines. 

To generate a representative LBG, we have started with the 
spectral energy distribution from Reddy et al. (2012). We add 
lines to match the observed EWs from Shapley et al. (2003); 
the objects shown in Fig. 21 in the bottom panel are represen- 
tative of objects with high Ly a EW in their sample. From 
these simulations it is clear that we will be able to detect both 
the Ly a and interstellar absorption features in the continua 
down to our detection limits. Finally the QSO spectra are 
generated from the composite spectrum from Vanden Berk et 
al. (2001). 

We also generate a picket-fence of emission lines across 
the entire waveband. We generate [O II] models as a pair of 
Gaussians at A, A 3726.2, 3728.9, with a ratio of 2 : 1. The 
pair of lines are assumed to have the same intrinsic width of 
°gas = 100 and 200 km s _1 for two different simulation. Ly 
a is modeled with a red-asymmetric profile. 

4.2.1. Redshifts 

Our greatest challenge comes in measuring absorption-line 
redshifts for z rj 2 galaxies at our flux limit. In practice, we 
expect that many of our targets will have some [O II] emis- 
sion, since most red galaxies have low EW emission lines 
(e.g., Ho et al. 1997; Yan et al. 2006). Nevertheless, there will 
certainly be passive red galaxies for which the 4000A break 
will be the primary source of redshift information. It will be 
measured in the NIR arm. The spectral simulations are very 
promising. It has been shown that reliable redshift measure- 
ments are possible if one achieves a signal-to-noise ratio of 10 
per pixel in the continuum around the 4000 A break (Tonry & 
Davis 1979). Thus, we expect that if we can characterize and 
subtract the sky as assumed here, we will be able to measure 
redshifts confidently at our continuum limit of Jab = 23.4 
mag. Furthermore, as shown by Chen et al. (2012), it is pos- 
sible to derive reliable stellar masses with 0.2 dex accuracy in 
spectra with S/N per pixel as low as 3, provided the 4000 A 
break and HS features are covered in the bandpass. 

4.2.2. Fluxes and Line-widths 
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FIG. 20. — Simulated spectra (see text) demonstrating the power of PFS. The spectra are plotted in units of "nanomaggies", a unit of flux used by the SDSS, 
defined such that one maggie is the flux from an AB= object, 3631 Janskys. Thus 1 nanomaggie is equivalent to an apparent AB magnitude of 22.5, or 
3631 X 10~ 9 Jy or 3.631 X 10~ 29 erg s — 1 cm -2 Hz -1 . The spectra have been binned to a resolution of R = 400, 300, 300 in the blue, red, and NIR anus 
respectively. Top Panel: This is a typical star-forming galaxy at z = 1.14 close to our detection limit, showing that we will have significant detections of [O II], 
H/3, and [O III], which we can use to measure gas-phase metallicity (Fig. 19). We show a three hour integration. These spectra are based on stellar population 
synthesis models of galaxies in the NEWFIRM Medium Band Survey. Bottom: A star-forming (left) and three passive (right) galaxies at z — 1.8 and magnitudes 
as indicated as seen in a three-hour integration with PFS. The spectra were also generated as model fits to the NEWFIRM Medium Band photometric data, which 
provides some constraint on emission-line EW. We show that at our magnitude limit we will have a significant detection of either [O II] or the 4000 A break for 
measuring redshifts. 



Many of our science cases depend on our ability to mea- 
sure not only the line centroid (the redshift) but also the line 
widths and fluxes. For [O II] lines with a combined flux of 
F = 10 erg s" 1 cm" 2 , we find that we can recover the 
total flux within < 20%, and the line-width to within < 40%, 
for all lines with S/N greater than 10 (over most of the band). 
We recall that this flux is typical of [O II] emitters at z « 2 
(Takahashi et al. 2007). We repeat the experiment for Ly a 
emission lines with F = 0.8 x 10~ 17 erg s" 1 cm" 2 , which is 
typical of LAEs at 5 < z < 7. Since we will be measuring Ly 
a in the blue and red arms, we can do a bit better, recovering 
fluxes and linewidths at the 10% level in the blue arm. In the 
red arm we can only recover linewidths at the 20% level (since 
the model is skewed, we fit a skewed Gaussian here). How- 
ever, these limits are for individual objects, while our primary 
goal is to stack at least 500 spectra in a given patch. Thus, 
we will do a factor of ~ 20 better at measuring the compos- 
ite Ly a profile of high redshift emitters. Finally, we attempt 
linewidth measurements for our QSO spectra. Even in half- 
hour exposures at our magnitude limit, we expect to measure 
the CIV linewidth with 10% precision. Our goal of measur- 
ing spectral shapes and estimating SMBH masses from these 
spectra will not be compromised by the S/N of the spectra. 
Obviously, there are well-known uncertainties in calculating 
reliable SMBH masses using broad UV lines, but we will be 
in a good position to perform internal cross-checks on various 
broad emission lines using simultaneous observations (e.g., 
Shen et al. 2008). 

4.3. Survey Design and Target Selection 



4.3. 1 . The Survey Fields 

Over the coming years, HSC will perform an ambitious 
imaging survey in grizy. The survey will have three com- 
ponents with different depths. That of most interest to us is 
the HSC Deep layer, which will reach a 5 er depth of i < 27.2 
mag. There is also narrow-band imaging planned that will be 
of crucial importance in selecting Ly a emitters at high red- 
shift. Of the more than 20 deg 2 covered by HSC Deep, we are 
most interested in the 16 deg 2 that already contain J— band 
data to at least J = 23.4 mag from the UKIDSS DXS survey. 

These fields (XMM-LSS, ELAIS-N1, and the COSMOS 
field) have been selected to maximize overlap with multi- 
wavelength data sets, including X-ray data to study accretion 
activity and mid-infrared data to probe obscured star forma- 
tion and accretion. In particular, we are excited about the 
prospect to use Spitzer and Herschel data for measuring com- 
plementary star formation rates to those derived from emis- 
sion lines or continuum fitting. For instance, we have good 
overlap with the Hermes fields (HerMES Collaboration et al. 
2012). Likewise, our survey fields will provide natural targets 
for ALMA and/or CCAT 

4.3.2. Sun'ey Design: The 1 < z < 2 Survey 

The survey designs are different for the continuum-selected 
galaxy sample at 1 < z < 2 and the study of high-redshift 
LBGs and LAEs at 2 < z < 7. For the 1 < z < 2 survey, the 
critical factors are: 

1 . We want to study the progenitors of local ~ L* galax- 
ies. Based on our spectral simulations, we can ex- 
pect to reach a S/N > 10 per resolution element (for 
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FIG. 21. — Top Panel: Similar to the previous figure, but for two QSO spectra assuming 30 minutes for the exposure time. The left-hand panel shows that for 
even a moderate-luminosity object at z = 7, a short exposure yields an unambiguous QSO with adequate signal/noise for measuring continuum and spectra line 
widths. Bottom: A star-forming galaxy selected as a {/—band (left) and g— band (right) dropout as observed in a six-hour integration with PFS. The spectral 



energy distribution is based on a spectral synthesis model from Reddy et al. (2012). We expect to detect Ly a with / 
level in 6 hrs of integration. Note we can detect absorption lines in individual spectra at J = 24.5 AB mag. 
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Fig. 22. — Le/h Illustrative models of evolving galaxies with different evolutionary histories based on Maraston (2005) stellar population 
models. All models have re _T star formation histories, and are normalized to be L* galaxies today. In red is a red galaxy progenitor with a 
final mass of M* = 6 x 10 10 Mq. In green is a more typical star-forming galaxy with a final mass of M* — 3.4 x 10 10 MQwith an extended 
history, while the blue line is a more extreme star forming system with a final mass of M* — 2.3 x 10 10 Mq. For reference, in black is 
the mean apparent J— band magnitude of a 3 x 10 Mq galaxy at each epoch from the A'— band selected NEWFIRM medium-band survey. 
Right: Expected stellar mass distribution for our 1 < z < 2 survey, based on the NEWFIRM Medium Band Survey (Whitaker et al. 201 1). We 
apply a J— band magnitude limit of 23.4 mag, but note that the total number of points in the figure correspond to a 0.2 deg 2 survey; we will 
observe 125 times as many objects. Without a redshift selection, a large fraction of the galaxies in a magnitude-limited sample have z < 1. The 
banding in the figure is due to large scale structure; our much larger volume will average over this cosmic variance. 



R ~ 300) in 3 hr exposures at our limiting magnitude 
of Jab = 23.4 mag. With this magnitude limit, we 
expect to probe masses as low as 10 10 M© at z rj 2 
based on the experience of NEWFIRM Medium-Band 



redshift survey (Fig. 22). Note that we are only com- 
plete for stellar masses > 10 11 Mq, since below this 
limit we miss the red, passively evolving galaxies. We 
have investigated a selection based on Spitzer/IRAC de- 
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TABLE 7 
Impact of Cosmic Variance 



TABLE 8 
All Targets 



Log M * 

9.1 

9.6 

10.1 

10.6 

11.1 

11.6 



1< z < 1.4 

0.02 

0.02 

0.02 

0.02 

0.03 

0.04 



1.4 < 2 < 1.8 

0.02 

0.02 

0.02 

0.03 

0.03 

0.05 



1.8 < 2 < 2.2 

0.02 

0.03 

0.03 

0.03 

0.04 

0.06 



2.2 < 2 < 2.1 

0.03 

0.03 

0.03 

0.04 

0.05 

0.07 



Mag. 



Selection N 



Int. Time Area 



NOTE. — The logarithmic error bars on the number density in each logarithmic 
stellar mass bin (log Mq) in the 16 deg 2 survey. Error bars are shown for four 
redshift ranges: 1 < z < 1.4, 1.4 < z < 1.8, 1.8 < z < 2.2, 2.2 < z < 2.6. 
The error bars are derived using the method of Moster et al. (2010). 

tection. However, we do not have deep enough IRAC 
imaging over all HSC Deep fields. Furthermore, we 
would have a large tail of galaxies with Jab > 23.4 
mag, for which we would be unable to measure a red- 
shift. Finally, selecting in the rest-frame optical ensures 
that we will observe both red and blue galaxies (Fig. 
22). Thus, we have chosen a rest-frame optical (ob- 
served J— band) selection as a compromise between 
UV-selection on the one hand, and stellar-mass selec- 
tion via Spitzer on the other. Further color selection 
will be needed to remove low redshift contaminants, as 
described in detail in point 4 below and §4.5.3. 

2. In order to observe the same population at low redshift, 
we also plan to observe all objects to a limiting magni- 
tude of Jab = 21 mag. Since we will spend only 20 
min per galaxy at these bright limits, and since we can 
rapidly reconfigure the fiber placement, we get these 
6000 galaxies per deg 2 essentially for free, and they 
provide a low-redshift comparison sample (e.g., for 
comparing Ha and [O ll]-derived star formation rates). 

3. The area should be large enough to minimize cosmic 
variance. In 16 deg 2 we cover a volume that is a few 
tenths that of the SDSS in a redshift bin at z « 2, which 
should sample a fair volume of the universe. To quan- 
tify the level of cosmic variance, in Table 7 we show 
the expected error bars on the number density of galax- 
ies in mass and redshift bins derived using the method 
of Moster et al. (2010). The nominal survey area of 16 
deg 2 will avoid limitations of cosmic variance. 

4. We want a sample selection that is as unbiased as pos- 
sible with regard to galaxy color, stellar age, line emis- 
sion, etc. We have considered various selection tech- 
niques, including photometric redshift selection, IRAC 
color, and optical-NIR color selection. None is perfect, 
but the pros and cons of each are outlined below. 

4.3.3. Target Selection: The 1 < z < 2 Survey 

A major strength of the SDSS galaxy survey was the simple 
magnitude-limited selection with an easy-to-model selection 
function. We would like to replicate that ideal to the extent 
possible. Recall that we must integrate for 3 hr per source at 
our magnitude limit of Jab = 23.4 mag to obtain scientifi- 
cally useful spectra. There are 12, 000 galaxies per deg 2 with 
23 < J < 23.4. Even with the multiplexing capabilities of 
PFS, a purely magnitude-limited survey over 16 deg 2 would 
take nearly 200 nights. Furthermore, > 50% of our targets 
would be faint galaxies at z < 1. As we select in redder bands, 
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NOTE. — Col. (1): J— band magnitude (AB) or narrow-band magni- 
tude. Col. (2): Redshift range. Col. (3): Selection method. Col. (4): 
Number of galaxies per deg 2 . For the continuum-selected sources, we have 
adopted the mean of the number densities from the COSMOS and NEW- 
FIRM catalogs. The LBGs are grouped in redshift with Sz R2 1. Col. (5): 
Integration time (hrs). Col. (6): Survey area (square degrees). 

the fraction of z > 1 targets increases. For instance, only 40% 
of 23 < J < 23.4 objects are low redshift contaminants (see 
also Table 8). We have chosen a J— band magnitude limit 
for this reason and because our spectra cover the J— band. 
Over the 16 deg 2 area, deep NIR imaging down to 23.4 mag 
is already available (in the fields XMM-LSS, ELAIS-N1, and 
the COSMOS field). Even with a J— band selection we will 
need to use additional color information to isolate the redshift 
range of interest. We will perform a purely magnitude-limited 
survey over two PFS pointings to be observed early in the sur- 
vey, which will allow us to calibrate our selection methods. A 
similar strategy was adopted by DEEP2 (Davis et al. 2003). 

We have explored various color-selection techniques to 
maximize the time spent in the redshift range of interest. We 
outline the strengths and weaknesses of three potential color 
selections here: 

1 . Photometric redshifts: As a straw-man, using photo- 
metric redshift is a very appealing strategy, as it utilizes 
all of the exquisite HSC imaging. We are testing the ca- 
pability of HSC Deep to measure reliable photometric 
redshifts using simulations based on COSMOS data. In 
the simulations, the completeness, or the fraction of ob- 
jects with z spec > 1 that correctly have z p hot > 1, is 
93%. The contamination, which is the fraction of ob- 
jects at z spec < 1 that have z p hot > 1, is 5%. Further- 
more, our ability to measure reliable photometric red- 
shifts should improve dramatically in the coming years, 
as we gain access to real HSC data, and take advantage 
of ongoing redshift surveys in our field (e.g., VIPERS). 
A major concern with using the photometric redshift is 
systematic error (for instance, in the template set). We 
will be measuring photometric redshifts for a popula- 
tion without good spectroscopic coverage, and there- 
fore there is room for significant spectral differences 
between the templates and the targets. Much more work 
is needed to determine whether photometric redshift se- 
lection can provide an unbiased sample selection. 

2. IRAC CH1-CH2 Selection: There is a peak in galaxy 
spectral energy distributions at A = 1.6 fim, red- 
ward of which is the Raleigh-Jeans tail. As this fea- 




Fig. 23. — Left: Differential redshift distributions for an i-band (AB) selection with the half-magnitude bins indicated at the top. These numbers 
indicate the lower magnitude of each bin. The number counts come from the 30-band photometric redshifts from COSMOS (Ilbert et al. 2010) 
over a 2 deg 2 area; the structures we see are due to cosmic variance. Right: Same, but for J— band (AB) selection. 



ture moves through the IRAC bandpasses (particularly 
CHI = 3.6 //mm and CH2= 4.8 pm), the CH1-CH2 
color goes through a distinct minimum at z « 0.8 
and then rises monotonically towards higher redshift 
(Fig. 24 left). Thus, a simple cut of CH1-CH2 > -0.1 
effectively selects z > 1.3 galaxies (e.g., Papovich 
2008). Using the NEWFIRM Medium Band Survey, 
we find that the completeness, or the fraction of ob- 
jects with z > 1 that are selected by this method, is 
64%, while the contamination, or fraction of selected 
objects that have z < 1, is 30%. The numbers are not 
as impressive as the photometric redshifts, but the se- 
lection is easily modeled and relatively free of system- 
atics. The problem is that we would have to obtain suf- 
ficiently deep Spitzer data (~ 2/iJy) over the entire 16 
deg 2 field to use this method. 

3. griJ Selection: Inspired by the DEEP2 color selection, 
we have designed an optic al/NIR color selection that 
is effective at removing z < 0.8 interlopers. Note that 
K— band imaging would be better for the purpose of ex- 
cluding low-redshift sources, but is not available. The 
color region we select is defined by i — J > 1.15 or 
g-r < (0.5 + 0.359[i - J + 0.8]) and Jab < 23.4 
mag (Fig. 24 right). Again based on the NEWFIRM 
Medium Band Survey, the completeness (fraction of 
galaxies with z > 1 that we select) in this case is 93%. 
The contamination (galaxies in the final sample with 
z < 1) is 42%. However, if we look at galaxies in 
the final sample with z < 0.8, the contamination is 
only 24%. Therefore, this simple color selection, with 
all bands to be observed as part of HSC or in hand, 
is very promising. With a redshift limit of z > 0.8, 
the total number of objects per deg 2 would increase by 
37% from the straw-man survey presented here. The 
numbers come from the NEWFIRM catalog but change 
little when we use the COSMOS photometric redshift 
catalog instead. In what follows, the numbers assume 
a perfect cut at z = 1, but the appropriate correction 
would only be an increase of 30% for the number den- 



sity of continuum-selected galaxies. 

We can further mitigate contamination through real- 
time redshift determinations. The main contaminant in a 
continuum-limited sample is moderate-redshift, intrinsically 
faint galaxies that are likely strongly star forming. We should 
be able to identify their redshifts with a single 20 min ex- 
posure, if we can carry out real-time analysis of the spectra. 
Thus, we can move to another target in < 1 hr, given the fact 
that the fibers are individually deployable. In addition, we 
will sharpen our target selection with our magnitude-limited 
component, as well as with the redshift information from the 
PFS cosmology survey and the VIPERS spectroscopy, which 
overlaps with our fields. 

If we make a perfect cut at z p hot ~ 1 and Jab < 23.4 
mag, then we expect 22000 galaxies per sq. deg between 
1 < z < 2. These numbers, based on small angle sur- 
veys, are still subject to cosmic variance at some level. From 
comparisons between NEWFIRM- AEGIS and NEWFIRM- 
COSMOS, large-scale structure could change our numbers by 
as much as a factor of two. In Table 8 we show the number 
of targets as a function of magnitude, number density, and ex- 
posure time, using the average number densities between the 
NEWFIRM and COSMOS catalogs. We plan for a total of 
68 nights for this component of the survey, assuming 3 hrs 
per galaxy in the faintest half-magnitude bin. In Fig. 22, we 
show that we will be able to observe galaxies with masses 
M >~ 10 10 Mq at our redshift limit, allowing us to study 
the progenitors of L* galaxies today. Thanks to the J— band 
selection, we will be sensitive to the red sequence even at 
z w 2, as shown in Fig. 25. The PFS survey is compared 
in redshift, depth, and area to other existing and upcoming 
surveys in Fig. 18. Our survey is designed to explore a poorly 
explored epoch in cosmic history, at the time when star for- 
mation was at its peak and galaxies began to resemble those 
we see today. 

Finally, we also propose to perform a purely magnitude- 
limited survey in 2 PFS pointings. With two independent 
pointings, we will sample different large-scale structures, and 
given that we expect ~ 10 clusters with M > 10 14 Mq with 
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Fig. 24.— Left: The dependence of Spitzer/IRAC CH1-CH2 (3.6 - 4.8^m) on redshift. Galaxies are from the NEWFIRM Medium Band 
Survey, and IRAC data have 1 fiJy depth. The structure is due to the 1.6/jm peak in galaxy spectra moving through the CHI band. Right: The 
gri J color selection, where galaxies are from the COSMOS survey. We select galaxies with i — J > 1.15 org — r < (0.5 + 0.359[z — J + 0.8]); 
the selection is highly complete (~ 90%) while the contamination is only ~ 20% for z < 0.8. 
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Fig. 25. — Rest-frame color-color plots of the NEWFIRM Medium Band Survey catalog, in the redshift bins indicated and with the J— band 
cut of 23.4 mag. In red are the color cuts used to define the red sequence by Williams et al. (2009). It is clear that we will sample the red 
sequence all the way to z ~ 2. 



z > 1, we do not expect to be cosmic variance limited. Ide- 
ally, this component would be completed early in the survey. 
This small (~ 10 night) investment, including observations of 
an additional 24k galaxies per deg 2 over the nominal color- 
selected survey, is essential for the following reasons: 

1 . It will allow us to tune and understand our color selec- 
tion in detail in preparation for the main survey. 

2. Our color cuts will certainly exclude certain rare popu- 
lations with extreme colors, such as luminous obscured 
quasars (Zakamska et al. 2003), massive post-starburst 
galaxies (Tremonti et al. 2007), and other strong line- 
emitting systems. Within a pure magnitude-limited sur- 
vey we can assess the missing populations, study their 
properties, and decide whether to design special color 
selection to target them. 

3. The bias of the galaxies used for the cosmology sur- 
vey may be a strong function of redshift, as it involves 



a color cut which may select astrophysically distinct 
populations of galaxies at different redshifts. The bias, 
and the relationship of the selected galaxies to the full 
galaxy population, needs to be understood in order to 
use the sample for cosmological inferences, especially 
when combining clustering measurements at different 
redshifts. (Note that the BAO scale itself is actually 
fairly insensitive to these details, as the bias is scale- 
independent on these very large scales). 

4.3.4. Survey Design and Target Selection: Galaxies at the 
Earliest Epoch 

For the 2 < z < 7 LBGs and LAEs the deciding factors in 
survey design are: 

• We require a wide enough area to study bright and rare 
Lyman break galaxies. Existing spectroscopic surveys 
are too narrow to study these populations robustly. Over 
our 16 deg 2 survey we expect to find 6400 U— band 
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dropouts brighter than i < 24 mag. 

• The area must be large enough to span multiple ionized 
bubbles in the IGM at z w 7, since the size of ionized 
bubbles are as large as ~ 1 degree at the end of the 
reionization epoch. 

We will identify LAEs and LBGs at 2 < z < 7 down to 
NB sa 25 mag and zab = 24 mag. We expect to spend 3-6 hrs. 
per target (Table 8) in the same ~ 16 deg 2 field. The LAEs 
will be selected from the HSC Deep and Ultradeep narrow- 
band imaging, while the LBGs will be selected as dropouts 
(Fig. 26). The narrow-band filters in the Ultradeep survey 
are NB387\z = 2.2, NB526;z = 3.3, NB7l7;z = 4.9, 
NB816; z = 5.7, NB921; z = 6.6, and NB101; z = 7.3. In 
the Deep survey, we will have imaging in NB387, NB816, 
and ^5921. We cannot possibly survey all LAE/LBG galax- 
ies over this area and to this depth; there are far too many 
objects. We will therefore target the following populations: 

1. At intermediate redshift (2.6 < z < 3.4) we will target 
the bright end of the LBG luminosity function down 
to Jab = 24 mag (Fig. 26) over the full 16 deg 2 
survey area. We will use a U— band dropout selec- 
tion, and the existing depth of U = 26 AB mag in 
the HSC Deep fields from the CFHT archive allows 
us to select dropouts down to i = 24.5 mag. We will 
also select riz— band dropouts with i < 24.4 mag and 
5 < z < 7. Such bright galaxies are relatively rare 
(~ 430 per deg 2 ) and so their intrinsic properties are 
not well-studied. These bright systems serve as a good 
tracer of large-scale structure; we will investigate their 
clustering properties, and their cross-correlation with 
quasars. In addition, the intrinsic evolution in Ly a lu- 
minosity from 2.5 < z < 4 will serve as an important 
benchmark for evaluating the impact of reionization on 
the Ly a luminosity function at z > 5. Taking the total 
number of targets multiplied by the exposure time and 
dividing by the 1700 science fibers per PFS pointing, 
we require 12 hours to execute this program. However, 
note that the targets are spread over the entire survey 
area and thus will not be observed in a continuous 12 
hr sequence. All exposure times quoted in this section 
have been calculated in this manner, and are only meant 
to demonstrate the fractional allocation to each survey 
goal. 

2. We will also target LBGs over the same redshift range 
to a slightly fainter limiting magnitude of zab = 24.5 
mag over the HSC Ultradeep area (3.5 deg 2 ). Our pri- 
mary goal is to use these targets as backlights to study 
the halo gas of intervening galaxies. We will stack the 
background systems as a function of impact parameter 
and study the kinematics of gas in the outer regions of 
these halos. We will probe galaxy halos at z « 1.3 us- 
ing the continuum selected sample at z ps 2, and halos 
at z pa 2 using background LBGs at z sa 3. The fore- 
ground sources will be part of the continuum-selected 
survey above, and thus will have measured redshifts. 
The number density of these sources is 18k per deg 2 , 
exposure times are 3 hr per galaxy, and we observe 
them over 3.5 deg 2 so we are allocating 19 nights to 
this program. 

3. Finally, we will target all LAEs with z = 5.7,6.6 to 
NB = 25 mag over the full 16 deg 2 area. We will 



use the evolution in line-shape, luminosity function, 
and clustering of this population as a probe of reioniza- 
tion. The wide area will allow us to study the patchiness 
of reionization that is expected to manifest on degree 
scales (e.g., McQuinn et al. 2007). These galaxies are 
very faint: we will devote 6 hrs per galaxy. Their num- 
ber density is 500 per deg 2 for both redshifts combined, 
so we need 30 hours for this component of the survey. 
Crucially, we will also study the intrinsic line shapes of 
LAEs at z = 2.2, for which we can obtain a systemic 
redshift from the [O II] line. It is important to quantify 
the intrinsic variability in Ly a line shape to distinguish 
from changes caused by neutral absorption at the epoch 
of reionization. We will observe these down to a limit- 
ing Ly a flux of 2 x 10 -16 erg s _1 cm -2 , corresponding 
to an [O II] flux of — 8 x 10~ 18 erg s _1 cm~ 2 . The ex- 
pected space density is ~ 50 per deg 2 but we will need 
to target ~ 100 per deg 2 for this success rate. We will 
have a final sample of 800 objects, and this will take six 
hours. 

4. In addition to these primary goals, we also have a few 
galaxy populations that may be targeted as filler to im- 
prove survey efficiency. Since each fiber has limited 
reach, we may be unable to find a suitable target for 
a given fiber at a given time. Auxiliary high redshift 
targets that we are considering include Ly a blobs (we 
expect ~ 150 in 16 deg 2 over the entire redshift range), 
LAEs at 3 < z < 5, and fainter LBGs at 2 < z < 5. 
In the latter cases there are many thousands of targets 
per deg 2 , but we will prioritize those in the vicinity of 
either quasar candidates or Ly a blobs. 

4.3.5. Target Selection: Quasars down the luminosity function 

Like the LBG samples described above, we will pursue a 
two-tiered quasar selection. However, since the quasar popu- 
lation is so rare, we will perform a shallow survey over 1400 
deg 2 and then a deeper survey over the 16 deg 2 HSC Deep 
area. Our target strategy is summarized in Table 8. 

The PFS-AGN wide survey is based on the wide layer of 
the HSC multi-band imaging legacy survey (5-band optical 
images of g, r, i, z, and y, combined with available wide NIR 
surveys by VISTA and UKIRT over 1400 deg 2 ). We focus 
on the color-selected quasar candidates at 3 < z < 7 us- 
ing color cuts similar to those used by the SDSS (Richards 
et al. 2002). We will target a few tens of quasar candidates 
per PFS pointing to a magnitude limit of i < 24, to achieve 
the number-density of quasars shown in Table 8. The primary 
goals include identifying 6 < z < 7 candidates for deeper 
spectroscopy to study reionization, measuring the quasar lu- 
minosity and mass function over 3.5 < z < 7, and measuring 
the quasar auto- and cross-correlation functions and their de- 
pendence on redshift, mass, and luminosity. The total number 
of targets is small, and it will be integrated into the wide- 
field cosmology survey. As shown in §4, we can achieve these 
goals in 0.5 hr exposure times even at our magnitude limit. 

The PFS-AGN deep survey is based on the deep layer of 
the HSC legacy survey. In addition to color selection, we 
will also target variable objects and X-ray and mid-infrared 
selected objects. We will go a magnitude deeper than the 
PFS-AGN wide survey, and expect to target 100 active galax- 
ies per deg 2 . The goals include finding ~ 10 faint z sa 7 
quasars, studying X-ray and variability selected active galax- 
ies at 1 < z < 4, cross-correlating absorption systems in 
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the quasar spectra with the foreground galaxy population, and 
studying obscured populations in the galaxy survey that are 
identified via narrow emission lines. If we spend an hour per 
target for the majority of these targets, then it will take very 
few fibers from the main galaxy evolution survey. 

4.3.6. Summary of the Survey 

In Table 8 we show all the sub-samples that we plan 
to target, their expected number density per square degree, 
the required exposure times, and the proposed survey area. 
Throughout, when we quote the required number of nights, 
we are assuming 1700 science fibers per PFS pointing and 
zero overhead in changing targets. We quote the number 
of required clear nights (that is, no weather factor has been 
applied). We outline a two-tiered survey, with our primary 
galaxy survey covering 16 deg 2 and a deeper component cov- 
ering 3.5 deg 2 To summarize, we will study the following 
galaxy populations: 

• We will perform a continuum-selected survey of 1 < 
z < 2 galaxies with Jab < 23.4 mag over 16 deg 2 
and < z < 1 galaxies with Jab < 21 mag. We 
will target ~ 28, 000 galaxies per deg 2 The primary 
survey will be color-selected, using a combination of 
HSC and IRAC imaging. We will also perform a purely 
magnitude-limited survey over two PFS pointings early 
in the survey to calibrate our selection strategy. 

• At 2 < z < 7, we trace the galaxy population using 
bright dropouts down to a limiting magnitude of i > 
24 mag over the full 16 deg 2 . We expect to target 430 
galaxies per deg 2 . We will survey fainter LBGs to i > 
24.5 mag over 3.5 deg 2 to study intervening absorption 
systems. 

• We focus on LAEs at z = 5.5, 6.6, 7.3, in order to use 
the Ly a emission as a probe of reionization. Addition- 
ally we will study LAEs at z ~ 2.3, where we can use 
the [O II] emission line to measure the systemic red- 
shift. 

Adopting these survey areas and target densities, the survey 
requires 105 clear nights. The survey parameters are summa- 
rized in Table 9. 

4.4. Requirements 

We are assuming the nominal sensitivities and resolution 
of PFS. The feasibility of our survey relies most heavily on 
the sensitivity of PFS in the red and NIR arms, from which 
we will derive redshifts and key diagnostic information. Our 
most important assumptions are: 

1. We assume 1700 fibers per PFS pointing. We also as- 
sume that it is possible to reposition fibers rapidly, so 
that we move off of bright targets in a short time. This 
large multiplexing capability is crucial for our science 
case, otherwise we cannot achieve the large areas and 
high sampling rates that make our survey unique. 

2. We assume that we can measure continuum-based red- 
shifts at our magnitude limit of Jab = 23.4 mag. 
Based on our spectral simulations, we find that in a 3 hr 
integration we will have at least the required S/N of 10 
per pixel in the NIR arm to measure continuum-based 
redshifts (Tonry & Davis 1979), assuming R ~ 300. 



We suggest real-time spectral analysis to identify low 
redshift interlopers after the initial 20 min integration. 

3. We expect to reach a 5 a emission-line limit of 8 x 
10 -18 erg s _1 cm -1 for > 50% of the band. This limit 
is reached in the NIR arm in 3 hrs of integration on the 
[O II] line. We are assuming intrinsic line widths in the 
galaxies of cr gas w 100 kms -1 . The S/N will degrade 
by ~ 30% for a galaxy with twice the dispersion, as the 
light is spread over roughly twice as many pixels. This 
gets us to J[oii] ~ 8 x 10 41 erg s~ 4 at z ps 2, which 
is below the [O II] luminosity of a typical L* galaxy at 
this epoch (e.g., Takahashi et al. 2007). 

4. We find that we can detect Ly a with 7 x 10 -18 erg 
s~ 4 cm~ 2 in the blue arm with R w 1300 at the 5cr 
level in 6 hrs of integration. Again we assume an in- 
trinsic gas dispersion of 100 kms - , which is typical of 
measured LAEs (Ouchi et al. 2010). We are thus de- 
tecting sources with L > 1.5 x 10 42 erg s -1 , which 
is a factor of five fainter than the typical LAE at these 
redshifts (Ciardullo et al. 2012). In the red arm, our 
limits are slightly better. In 6 hrs our 5 a detection limit 
is 5 x 10~ 18 erg s _1 cm~ 2 for > 50% of the band. 
Here our limits of 8.5 x 10 42 erg s~ 4 at z = 5 and 
2 x 10 43 erg s _1 at z = 6 are around or slightly ex- 
ceed L* Lya at these epochs (Ouchi et al. 2010). Our 
limits are comparable to those achieved by Ouchi et al. 
(2010) for an individual object, and we will increase the 
sample size by two orders of magnitude. 

5. In order to achieve the above continuum and emission- 
line limits above, we place strong requirements on our 
software, specifically in regard to sky subtraction in the 
red and NIR arms. We currently assume that we are 
able to determine the sky level at a given time to 0.5%. 
This is easily supported by the Poisson noise in the 
lines themselves. We further assume that we achieve 
a wavelength solution that is good to 0.01 pixel. Fi- 
nally, we assume that no more than 1% of the light hit- 
ting the detector is distributed across the entire chip as 
a smooth background. We therefore assume that we un- 
derstand systematics in sky subtraction at this level, and 
that there are no other systematic errors. 

5. SCIENCE REQUIREMENTS 

We now summarize the science requirements arising from 
the earlier sections. These requirements are necessary for de- 
riving lower level instrument requirements and in understand- 
ing the flow of requirements from the science to detailed en- 
gineering. Establishing the science requirements flow is like- 
wise essential in order to track the science impact of potential 
changes to the instrument configuration in the future. The top 
level science requirements and observing requirements are in- 
tended to be directly verifiable items in the project. For con- 
venience, we list the overall drivers, item by item, both in 
tabular form and in more detail below. 

The top-level science requirements for the PFS surveys de- 
scribed in the previous chapters are summarized in Table 10 22 . 

We briefly discuss the individual entries in Table 10. 

22 Our instrumental requirements for the different scientific goals outlined 
in this document are somewhat different. In particular, the galaxy science, 
which requires the deepest exposures, have the strongest requirements on 
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TABLE 9 
Survey Parameters 



1 < z < 2 

Selection (optical) HSC Deep: i < 27.2 mag 

Selection (NIR) XMM-LSS, ELAIS-N1, and COSMOS: J AB < 23.4 mag 

Exposure time at flux limit 3 hrs 

Area 16 deg 2 

Number of nights 68 nights 
Hux-limited survey 

Selection (NIR) ' XMAf-LSS, ELAIS-N1, and COSMOS: J AB < 23 A mag 

Exposure time at flux limit 3 hrs 

Area 2 PFS pointings 

Number of nights 10 

High-redshift galaxies 

Selection (optical) HSC Deep: i < 27.2 mag 

Selection (narrow-band) HSC Deep: NB < 25 mag 

Exposure time at flux limit 3 hrs 

Area 16 deg 2 

Number of nights 27 

Total number of nights 105 



3.0 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 p-n-i 




Fig. 26. — Left: Color magnitude diagram of narrow-band excess color (z — NB921) vs. narrow-band magnitude (7V.B921) for selecting 
z = 6.6 LAEs in the SXDS field (Ouchi et al. 2010), which will be applied to the LAE selection of our PFS survey. Black dots present colors of 
all the detected objects, while black filled circles and triangles denote the z = 6.6 LAEs. Red and blue open symbols mark spectroscopically- 
identified objects in the redshift range of LAEs (z — 6.45 — 6.65) and interlopers, respectively. The green lines indicate 2a errors of the 
color of z' — NB921. Dashed and dotted lines represent the color cut for the narrow-band excess and the la limit of z' data, respectively. 
Right: Two color diagram of objects down to i' < 26.85 to select z ~ 4 LBGs in the Subaru Deep Field (Yoshida et al. 2006). The horizontal 
sequence along B — R = 0.62 shows objects which are fainter than the la magnitude in both B and R. The colored symbols show objects 
with spectroscopic redshifts, where cyan, green, violet, and red represent objects in the range z < 3.0, 3.0 < z < 3.7, 3.7 < z < 4.0, and 
4.0 < z < 4.5, respectively. The thick orange line indicates the boundary which are adopted for the selection of z ~ 4 LBGs. Our PFS survey 
will use a broad-band set of g'r'i' different from this BRi' bands for selecting z ~ 4 LBGs. However, the contamination rates of our z ~ 4 
LBG sample will be comparable, as demonstrated by another z ~ 4 LBG selection with GRI bands (Steidel et al. 1999). 



• Wavelength Coverage: The triple-armed spectrograph 
design offers complete and simultaneous spectral cov- 
erage from 380 through 1300 nm and clearly this rep- 
resents a unique feature of PFS. The driver for this re- 
quirement arises from all 3 surveys discussed in Chap- 

systematics such as stray light and sky subtraction. The overall requirements 
listed in Table 10 represent a compromise of sorts between the needs of the 
different scientific goals. As the instrument design becomes more solid, we 
will see to what extent the more stringent requirements of the galaxy evolu- 
tion science are possible. 



ters 2-4 but particularly the cosmology and galaxy sur- 
veys. The cosmology survey will enable the first BAO 
and RSD tests beyond z ~2 and the galaxy survey will 
track emission line properties of an unprecedented sam- 
ple of galaxies with no 'gaps' in redshift coverage. 

• Spectral Resolution: The resolutions of the three arms 
of each spectrograph are designed to ensure optimal ef- 
ficiency in the study of stars and galaxies over the wide 
wavelength range. The resolution must be higher in the 
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TABLE 10 
Top-Level Science Requirements 



Requirement 


Value 


Main Driver(s) 


Wavelength Coverage 


380- 1300nm 


All surveys 


Spectral Resolution 


Blue: R ~ 2500 
Red: R ~ 3000 
NIR: R ~ 4000 


Matched survey efficiency 

sky subtraction, stellar metallicities 

radial velocity precision 


Fiber Density 


> 2400 per 1.3 deg. diameter 
hexagonally-shaped field 


n a P g ~ a few @k = 0.1 h/Mpc 
requirement in BAO survey 


Throughput 


Average: >22% (blue) 
>22% (red) 

> 24% (NIR) 

Worst part of band: >20% (blue) 
>20% (red) 

> 18% (NIR) 


Matched survey efficiency 


Fiber Reconfig. Time 


£ 3min in total 


15min exposure in BAO survey 


Sky Subtraction Accuracy 


£j 1% of sky background 
per 4-pixel resolution element 


in all surveys 


Stray Light 


£j a few % of total sky brightness 
spread over detector 


in all surveys 


Read Out 


< 3 (Red) e~ rms per pixel 

< 4 (NIR) e" 


in all surveys 


Wavelength Calibration 


< O.lA 

(~l/20th of a resolution element) 


for velocity precision 


Fiber Diameter 


1.13" deg. diameter at the field center 
1.03" at the edge 


Optimal S/N for galaxy survey 


Pre-Imaging Data 


2 

gri to 25mag over 1400 deg 
griz+NB515 survey to V ~ 22.5mag 
grizy+NBs+J to ~ 27mag over 16 deg 2 


Color selection of ELGs for BAO survey 
Selection of stars in Milky Way/M3 1 
Color selection for galaxy survey 


Survey Area 


~1400 deg 2 for BAO 
~500 deg 2 for GA survey 
~ 16 deg 2 for galaxy survey 


Statistical requirements, cosmic 
variance, matching HSC data 



a Excludes atmosphere, central obscuration + WFC vignetting, and fiber aperture effect. 



red and near-infrared arm to ensure optimal sky sub- 
traction in the dense OH forest. The minimum effective 
resolution in the red arm is set by our requirement to 
measure accurate velocities and stellar metallicities for 
the Galactic surveys. 

Fiber Density: The fiber density is set by the field 
of view of the provided Subaru prime focus correc- 
tor and our requirement that we optimally sample the 
0.8 < z <2.4 galaxy population for the cosmology sur- 
vey at the wavenumber k where the BAO signal is to 
be recovered. Formally, this sets a requirement that 
the mean number density of target galaxies n g and the 
galaxy power spectrum P g (k) should satisfy n g P g ~1 
(or a few) at k ~0.1 h^ 1 Mpc. 

Throughput: The throughput requirement is driven by 
a combination of the S/N required for the various sur- 
veys and the available observing time, and is consistent 
with a preliminary assessment of the instrument per- 
formance. The values in the table exclude atmosphere, 
central obscuration, vignetting, and the fiber aperture 
effect. 

Fiber Reconfiguration Time: The reconfiguration time 
is set by the cosmology survey which is likely to have 



the minimum exposure time given it must survey a large 
volume by targeting strong line emitting galaxies. If the 
exposure time is ~15 minutes, a reconfiguration time 
(including slewing, acquisition, guiding, fiber position- 
ing and verification) should be completed in under 3 
minutes to minimize survey inefficiency. 

• Sky Subtraction Accuracy: This is set by the require- 
ment for high quality spectra at faint limits in all sur- 
veys. It is particularly driven by the need for accurate 
absorption line measures in the red and near-infrared 
spectral regions where the OH forest is dense. 

• Stray Light: As for sky subtraction, given the intensity 
of the OH night sky emission. 

• Detector Read Out: This is set by the signal/noise re- 
quirements of each survey, noting their likely photon 
and sky-noise components. 

• Wavelength Calibration: This is determined by the ra- 
dial velocity accuracy for the Galactic surveys. 

• Fiber Diameter: This is determined by the optimal sig- 
nal/noise requirements of the galaxy and cosmology 
survey, noting the distribution of line emission in these 
extended sources. 
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• P re-Imaging Data: This is determined by each individ- 
ual survey and the various color-criteria to be adopted 
in optimally selecting either galaxies in a particular red- 
shift range with certain star formation characteristics or 
stars of a given luminosity class in the Milky Way or 
M3 1 . The aim is to ensure that all the relevant data not 
currently in hand is taken by HSC. 

• Survey Area: This varies from one survey to another but 
is driven by statistical and/or cosmic variance require- 
ments and the availability of HSC imaging data. 

6. OUTSTANDING ISSUES 

Our primary goal has been to present the science case for 
PFS in the context of an anticipated 300-night Subaru Strate- 
gic Program (SSP) and, thereby, to define the scientific re- 
quirements for the instrument as summarized in the previous 
section. The successful review of this science case and its as- 
sociated technical document has led to an agreement to begin 
PFS construction with an anticipated first light in early 2017. 

Science planning will naturally continue through to the Pre- 
liminary Design Review (anticipated in early 2013) and here 
we briefly describe those activities in the context of instru- 
mental choices, the design of the survey and the broader im- 
pact of PFS beyond the initial SSP. 

It should be understood that PFS planning involves two 
rather different communities: the international partnership fo- 
cused scientifically via the PFS Survey Committee, and the 
Japanese astronomical community, whose interests are co- 
ordinated by the National Astronomical Observatory of Japan 
and its appointed national committees. 

In the following we briefly introduce some of the items 
which will form the basis of the project discussions prior to 
the PDR. 

6.1. Requirements Flowdown 

As discussed in Section 5, for a project with the size and 
complexity of PFS, formal science and instrument require- 
ments are necessary to guide the design and construction of 
the system, and allow the rapid assessment of design or per- 
formance changes during the course of the project. They will 
also guide the instrument test plan to be developed at a later 
stage in the project; each requirement will include the method 
for validating that it has been met. The full documentation 
of formal and tracked requirements will be completed by the 
PDR, and will define the required performance of the instru- 
ment system and the subsystems in order to achieve the sci- 
ence goals. 

The requirements can be broken into several levels, and we 
will use the following definitions. The Level 1 (LI) require- 
ments are the top level science requirements that comprise 
the reason for building the instrument. The Level 2 (L2) re- 
quirements are ones that are imposed on the entire instrument 
or multiple subsystems. Requirements on a single subsystem 
are Level 3 (L3) requirements. Level 4 (L4) requirements are 
the requirements created by the subsystem designers in order 
to implement the Level 3 requirements. 

Although the top level science requirements, L2, are sum- 
marized in Section 5, they will be refined and formalized by 
PDR (the LI requirements are described in detail in Sections 2 
- 4 for each science case). These focus on the science driven 
requirements and may not yet be complete. Specifically, they 
do not yet include all of the instrument level operational and 
interface requirements. The L3 and L4 requirements are in the 



process of being developed. Although there has been exten- 
sive discussion of these requirements, they are not yet ready 
for this documentation. 

The PFS project office will lead the effort to refine, for- 
malize, and document the L1-L4 requirements, in coordina- 
tion with the science and instrument working groups. This 
process will be complete for the L1-L3 requirements before 
PDR, while the definition of L4 requirements may evolve af- 
ter the PDR. The documentation will consist of a formal list of 
individual requirements with a description, rationale, and val- 
idation method listed for each requirement. The rationale will 
consist of a short paragraph explaining why the requirement is 
needed, any assumptions that were made, and what design ef- 
fort drove the requirement. It puts the requirement in context, 
and helps reduce ambiguities and misunderstandings. All of 
the instrument requirements are intended to be verified during 
the acceptance testing of the instrument, under a written test 
plan. 

6.2. Technical Design 

The outstanding item here is related to an intermediate dis- 
persion requirement which would significantly enhance the 
ability of PFS to measure chemical abundances for a wider 
range of elements. The Galactic Archeology working group 
has made a strong case for an additional spectroscopic mode 
with spectral resolution i?=5000 (and possibly higher) but 
with a limited wavelength coverage. This would enable PFS 
to unravel the temporal chemical evolution of each building 
block of Galactic halo and thick disk through their multi- 
element alpha abundances revealed in this wavelength range, 
namely [X/Fe] with X=Mg, Si, Ca, Ti in comparison with 
[Fe/H]. 

The PFS team has acknowledged the strong scientific case 
for this component which is also well-suited to the Sub- 
aru 8.2 meter aperture giving PFS a distinct advantage over 
massively-multiplexed spectrographs on 4 meter telescopes. 
Incorporating this requirement into the design was not possi- 
ble by the CoDR and we intend to evaluate and cost two pos- 
sible technical solutions prior to PDR as well as to determine 
their operational feasibility. 

• Reconfigure the current red arm spectrograph so that it 
is possible to switch from the default mode to one suit- 
able for intermediate dispersion studies. Such an 'inter- 
mediate dispersion mode' might be undertaken during 
a dedicated period of Galactic science or at the end of 
the cosmology/galaxy surveys. 

• Develop an independent spectrograph whose resolution 
is not restricted to lie within the boundaries of the cur- 
rent 3-arm design. Although a more expensive op- 
tion this may better suit a dedicated Galactic campaign 
given the earlier WFMOS studies made a strong case 
for a much higher resolution instrument. 

6.3. Planning the Subaru Strategic Program 

The three Working Groups have, thus far, largely worked 
independently in defining their science requirements and there 
has been little discussion of how to integrate the 3 surveys 
into a single coherent observing program. Ultimately the PFS 
partnership will have to take Sections 2-4 forward and write 
one or more Subaru Strategic Programs for consideration by 
the NAOJ Subaru Time Allocation Committee. 

Issues that the PFS team will address in discussions with 
NAOJ and the HSC team will include: 
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• Selecting fields within the HSC imaging survey 



• Science prioritization of the various survey components 



• Target sharing and exposure time optimization 



• Bright and dark time balance on the telescope 



The bright and dark time balance will depend critically on 
the question of how the Subaru Observatory intends to oper- 
ate PFS alongside HSC and the bright time use of PFS will, in 
turn, depend on whether there is a feasible intermediate dis- 
persion capability. 



6.4. Legacy Value 

The Subaru archive is formally a responsibility of NAOJ 
but the quality and delivery of the survey products is that of 
the current PFS science team. No negotiations have yet taken 
placed on this important issue, but we expect a full agreement 
on responsibilities, proprietary periods, international access 
etc by the time of the PDR. 



6.5. Synergy with Other Facilities 

PFS will have a significant impact beyond the surveys dis- 
cussed in this document in two ways. 

Firstly, the survey fields will deliver spectroscopic data of 
value to extant imaging facilities. During the pre-LSST era we 
can expect inquiries from those international teams for access 
to PFS data. Such access could be dealt with largely under 
the auspices of the Subaru archive discussed above, but there 
may be cases where new collaborations are valuable, either to 
the PFS team, the Japanese community or both. Such oppor- 
tunities will need to be considered and the roles and rights of 
the PFS team, NAOJ etc fully understood. 

A potentially more interesting opportunity is the likelihood 
of new collaborations with both current imaging survey tele- 
scopes and future facilities such as LSST, Euclid and perhaps 
ultimately TMT, beyond the scope of the Subaru Strategic 
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Program (SSP) introduced here. Formally, beyond the SSP, 
PFS will become the property of NAOJ but the current part- 
nership may have the opportunity to participate and contribute 
in future surveys conducted in conjunction with other facili- 
ties. 



7. SUMMARY 

Our intention via this document is to describe in as much 
detail as practical, the exciting scientific programs we have 
designed for the Subaru Prime Focus Spectrograph, as well as 
to present our first iteration for a Subaru Strategic Program 
(SSP) of ~300 nights that would realize these goals. The 
flow-down from these science requirements has been summa- 
rized in the context of the top-level instrument requirements 
which provide important constraints on the technical design 
of each component of PFS. 

PFS will likely be the first massively-multiplexed spectro- 
graph on a large aperture telescope to achieve first light. Its 
unique capabilities on an 8 meter platform at an excellent site 
promises exciting discoveries in 3 broad areas: the nature of 
dark energy through an ambitious survey of emission lines 
galaxies over an unprecedented redshift range 0.8 < z <2.4, 
the assembly history of the Milky Way and M31 through 
strategically-designed programs that exploit astrometric data 
from the Gaia mission, and galaxy evolution over a wide range 
in redshift (1< z < 7) exploiting PFS' unique wide wave- 
length coverage. 

Although our initial goal is to construct the instrument and 
conduct a refined version of the SSP presented in this article, 
we are strongly motivated to consider the longer term role of 
this survey instrument in the era of LSST, Euclid and TMT. 
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